NMR Quantum Computation 

Jonathan A. Jones 

Centre for Quantum Computation, Clarendon Laboratory, Parks Road, Oxford 0X1 3PU, UK and 

Oxford Centre for Molecular Sciences, New Chemistry Laboratory, South Parks Road, Oxford 0X1 3QT, UK 

(Received 31st July 2000) 



Contents 

Introduction 

Structure and scope 

Limits to computation 

Computational complexity 
Quantum complexity 

Atomic computation 

Computational circuits 
Keversiblc computation 
Reversible function evaluation 

Quantum computat ion 

Quantum parallelism 



Tomography 



Deutsch's algorithm 

Quantum gates 

Universality ot quantum gates 

Building JNMK quantum computers 

Qubits 



Logic gates 

Initialisa tion 

Readout 

Some two spin systems 

Scaling the system up 

Qubits and JNiVlK spin states 

Une qubit states 
'i'wo qubit states 



NMR logic gates 

One qubit gates 
Controlled- IN U'l' gates 
Uther two qubit gates 
Cates in larger spin systems 

Multi qubit logic gates 

Single transition selective pu lses 
Ceometric phase-shift gates 

Initialisation and NMR 

Spatial averaging 
Logical labeling 
'iemporal averaging 
I'he use ot "cat" states 

Readout 



Simple readout 



8 
8 
8 
9 
9 
10 

10 
10 
11 

12 
12 
12 
13 
13 
13 
14 
14 

14 
15 
15 
15 
16 

17 
17 



Practicalities 

Selective puhes 



Composite 1 



Abstract reterence frames 



Simple algo'ithms 



inunctions ar d phases 



Deutsch's algorithm 



IN MR implci icntations 



Extensions and simplihcations 



"G? 



over s qua 



antuni search 



IN MR implci icntations 



Lxtensions 



Quantum plienomena 



Lntaiio led st ates 



Quantum te. cportation 



Lrror corrcc ion 



NMR and e 



ntanglement 



Quantifying entanglement 



NMR and qi.antum mechanics 



Conclusions 



Acknowledg ements 



pseudo-entangled state 



References 



ulses 



I An explicitljT separably decomposition of a 



17 

18 
18 
18 
18 

19 
19 
19 
20 
21 
22 
23 
23 

24 

24 
24 
25 

26 
26 
27 

28 

28 



28 
29 



INTRODUCTION 



long before more formal publication. 



Quantum computation [0, ||, §] offers the prospect of 
revolutionising many areas of science by allowing us to 
solve problems well beyond the power of our current clas- 
sical computers [^. In particular a quantum computer 
would be superb for simulating the behaviour of other 
complex quantum mechanical systems ra, o. Although 
the theory of quantum computation has been studied for 
many years, and many important theoretical results have 
been obtained, early attempts to actually build even the 
smallest quantum computer proved extremely challeng- 
ing. 

In recent years there has been considerable interest in 
the use of NMR techniques fflto implement quantum 
computations ||, ^, |l^, ^ [l^, |13|, Q . It has proved sur- 
prisingly simple to build small NMR quantum comput- 
ers, and while such computers are themselves too small 
for any practical use, their mere existence has brought 
great excitement to a field largely deprived of experimen- 
tal achievements. 



Structure and scope 



LIMITS TO COMPUTATION 

Over the last forty years there has been astonishing 
progress in the power of computational devices using 
silicon based integrated circuits. This progress is sum- 
marised in a set of "laws" , usually ascribed to Moore, al- 
though some were developed by other people. Over this 
period the power of computational devices has roughly 
doubled every eighteen months, or increased ten fold ev- 
ery five years. Unfortunately, this extraordinary tech- 
nological progress may not continue for much longer, as 
this increase in computing power requires a correspond- 
ing decrease in the size of the transistors on the chip; this 
shrinking process cannot be continued indefinitely as the 
transistors will eventually be reduced to the atomic scale. 
At current rates of progress it is estimated jl^l that the 
ultimate limits of this approach will be reached by about 
2012, and any further progress in the power of our com- 
puters will require a radically different approach. One 
possible approach is to use the power offered by quan- 
tum computation. 



In this article I will describe how NMR techniques 
may be used to build simple quantum information pro- 
cessing devices, such as small quantum computers, and 
show how these techniques are related to more conven- 
tional NMR experiments. Many tricks and techniques 
well known from conventional NMR studies can be ap- 
plied to NMR quantum computation, and it is likely that 
some of these will be applied in any large scale quantum 
computer which is eventually built. Conversely, tech- 
niques from quantum computation could have applica- 
tions within NMR. 

It is impossible to explain how NMR may be used to 
implement quantum computations without first explain- 
ing what quantum computation is, and what it could 
achieve. This in turn requires a brief discussion of classi- 
cal reversible computation |]15| . In order to reduce these 
introductory sections to a reasonable length many tech- 
nical points have inevitably been skipped over. Wherever 
possible I will use traditional product operator notation 
pql to describe how NMR quantum computers are im- 
plemented; it will sometimes, however, be necessary to 
use more abstract quantum mechanical notations u^ to 
describe what these computers seek to achieve. 

Before the advent of NMR quantum computation, al- 
most all research in this field was performed within 
a small community of physicists and mathematicians. 
Some important results have never been published as 
conventional papers, but simply circulate as electronic 
preprints; copies of most of these can be obtained from 
the quant-ph archive on the LANL e-print server p8J. 
Similarly, many other papers appear as LANL e-prints 



Computational complexity 

Even if the problems described above are sidestepped 
in some way, there are still strong limits to the problems 
which we can solve with current computers. These limits 
are derived from the underlying theory of computation 
itself, and in particular from computational complexity 
theory [g^. Complexity theory considers the classifica- 
tion of mathematical problems according to how difficult 
they are to deal with, and seeks to divide them into those 
which are relatively easy (tractable) and those which are 
uncomfortably difficult (intractable). Note that complex- 
ity theory is not concerned with problems which we do 
not know how to solve, or which we know cannot be 
solved, but only with problems for which an algorithmic 
solution (tractable or intractable) is known. 

The classical theory of computation has remained 
largely unchanged for decades, with a central role played 
by the Church-Turing thesis. This asserts that all phys- 
ically reasonable models of computation are ultimately 
equivalent to one another, so that it is not really neces- 
sary to consider any particular computer when assessing 
the complexity of a problem; any reasonable model will 
do. In particular, the tractability or otherwise of a prob- 
lem is independent of the computational device used. As 
we will see, quantum computation challenges this thesis, 
as quantum computers appear to be fundamentally more 
powerful than their classical equivalents. 

To make further progress it is necessary to use a more 
precise measure of the complexity of an algorithm. The 



usual approach is to determine the computational re- 
sources (most commonly defined as the number of el- 
ementary computational operations) required to imple- 
ment it. Clearly this measure will depend on the exact 
nature of the computational resources available to our 
computer, and so this is not directly useful. A better ap- 
proach is to consider not a single isolated problem, but a 
family of closely related problems, and to determine how 
the resources required scale within this family. To take a 
simple example, adding two n digit numbers with pencil 
and paper requires n separate additions, together with n 
carries, and so the time required for addition scales lin- 
early with n. Similarly, multiplying two n digit numbers 
by long multiplication requires n^ multiplications and a 
similar number of additions. 

Mathematically, adding two n digit numbers is said 
to be 0{n) (that is, of order n). The exact number of 
steps required will depend on exactly how elementary 
computational steps are counted, but the total number 
of steps will always scale linearly with n. Similarly, long 
multiplication can always be performed in a number of 
steps which varies quadratically with n, and so is 0{n'^). 

The complexity of a problem, rather than an algo- 
rithm, may be defined as the complexity of the best 
known algorithm for solving the problem; thus addition 
is 0{n). It might seem that multiphcation is 0{'n?), but 
in fact long multiplication is not the best known algo- 
rithm; a better approach is known with complexity about 
O(nlogn) iQ. Strictly speaking one should also distin- 
guish between an upper bound on the number of steps 
known to be sufficient to solve a problem (indicated by 
O) , and a lower bound on the number of steps known to 
be required to solve a problem (indicated by fi); a prob- 
lem, such as addition of n digit numbers, which is both 
0{n) and Vt{n) is denoted as Q{n). 

Problems, such as addition and multiplication, whose 
complexity is at worst a polynomial function of n, are 
said to be easy, while problems whose complexity is worse 
than a polynomial function of n are said to be hard. For 
example, consider the problem of finding the prime fac- 
tors of an n digit composite number. The obvious way to 
do this is to simply try dividing the composite number 
by every number less than its square root; as there are 
approximately -s/TO" = 10"/^ such trial divisors, this al- 
gorithm has complexity 0(10"^^). Better algorithms for 
factoring are known, but they all have the same property 
of exponential complexity. 

For problems with polynomial complexity, especially 
those whose complexity is described by a low power of n, 
such as n or n^ , the effort required to solve a problem in- 
creases only slowly with n. Thus it should be possible to 
tackle such problems for a reasonable range of input sizes, 
and a modest increase in computer power should give a 
significant increase in this range. Exponential functions, 
however, rise extremely rapidly with n, and so it will only 
be possible to solve problems with exponential complex- 



ity for relatively small input sizes; furthermore a modest 
increase in computer power will result in only a tiny in- 
crease in the range of problems that can be solved. 

The apparent exponential complexity of factoring is a 
matter of some importance, as it underlies many cryp- 
tographic schemes in use today, notably those based on 
the Rivest-Shamir-Adleman (RSA) approach to public 
key cryptography, such as PGP (Pretty Good Privacy) 
pl[ . These schemes involve a public key, which anyone 
may use to encrypt a message, and a private key, which 
is required to decrypt such messages. The relationship 
between these two keys is equivalent to the relationship 
between the product of two large prime numbers and the 
two prime numbers themselves. The security of the sys- 
tem depends on the difficulty of determining the private 
key from a long public key, which itself depends on the 
complexity of factoring. By contrast, the computational 
complexity of encrypting and decrypting messages is only 
a polynomial function of the size of the keys. Thus a 
small increase in the amount of effort required to use the 
cryptographic scheme results in an enormous increase in 
its security. 



Quantum complexity 

Quantum computation offers the possibility of bypass- 
ing some of the limits apparently imposed by complexity 
theory. This is because a quantum computer could im- 
plement entirely new classes of algorithms which would 
allow currently intractable problems to be solved with 
ease. 

The first serious discussion of quantum computation 
was by Feynman, who analysed the difficulty of simu- 
lating quantum mechanical systems using classical com- 
puters [gl. This difhculty is well known and easy to un- 
derstand; it arises from the enormous freedom available 
to quantum systems. For example, a system of n cou- 
pled spin-i nuclei inhabits a Hilbert space of size 2", 
and so must be described by a vector with 2" compo- 
nents. Thus it appears that any classical algorithm to 
simulate the behaviour of n spin-^ particles must have 
complexity at least 0(2"); within NMR this corresponds 
to the well known computational difficulty of simulating 
a large strongly coupled spin system. This apparently 
inevitable exponential complexity makes the simulation 
of quantum mechanical systems an intractable problem. 

Despite this apparent complexity, however, coupled 
spin systems evolve in the "correct" manner. Thus, in 
some sense, such a spin system appears to have the ca- 
pacity to solve a problem which is intractable by conven- 
tional classical means. Clearly using a system to simulate 
itself is not a huge step forward, but Feynman suggested 
that it might also be possible to use one quantum me- 
chanical system to simulate another quite different sys- 
tem. Thus an easily controllable system might be used 



to simulate the behaviour of another less well behaved 
system, while a carefully chosen system might be usable 
as a general purpose quantum simulator ^, g] . 

Feynman's ideas appear to have been limited to the 
simulation of physical systems, and the ideas he de- 
scribed have more in common with analogue computers 
than with current digital computers. In 1985, however, 
Deutsch extended these ideas and described a quantum 
mechanical Turing machine, which could act as a gen- 
eral purpose computer |22 . Deutsch also described a 



(somewhat contrived) problem [M, E4l E5| which could 
be solved more efficiently on a quantum mechanical com- 
puter than on any classical computer, suggesting that it 
might be possible to use a quantum computer to solve 
otherwise intractable problems. 

Since that time several quantum algorithms pq | have 
been developed, the most notable of which is Shor's quan- 
tum factoring algorithm El ; this allows a composite num- 
ber to be factored with a complexity only slightly greater 
than 0{n^), thus rendering factoring tractable. As well 
as being of great mathematical interest, this algorithm 
has obvious practical significance, as it poses a threat 
to many current cryptographic schemes. The discovery 
of Shor's algorithm triggered an enormous increase in 
research directed at quantum computation and related 
areas of quantum information theory. 




FIG. 1: A circuit to compute the exclusive-OR (xor) func- 
tion, z — X XOR y, using AND, OR and NOT gates. Note that 
this circuit uses three implicit gates, two CLONE gates, shown 
as small circles, where wires split into two (to copy the in- 
put) and one SWAP gate (where the wires cross over). These 
implicit gates are fairly easy to implement in traditional elec- 
tronic computers, but can cause problems in other designs 
and so cannot simply be ignored. Some authors even con- 
sider the wires which interconnect gates as non-trivial gates 
in their own right. 
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ATOMIC COMPUTATION 



Before describing how quantum mechanical computers 
can be built, it is useful to consider how an atomic scale 
system, such as a group of coupled nuclei, could be used 
to implement classical computations uB, E^, p8[ . Discus- 
sions of this predate suggestions that quantum devices 
might have fundamentally different properties. 

The basic approach is very simple. Classical computa- 
tion [ p5[ is implemented with two state devices, usually 
called bits, and these two states can be mapped to the 
two eigenstates of a quantum mechanical two level sys- 
tem. The calculation then proceeds by manipulating the 
states of various bits such that the final state of some 
group of bits (the "output register" ) depends on the re- 
sult of the desired computation. As will be shown later, 
quantum computation is very similar, except that the 
bits are not confined to their eigenstates; this effectively 
allows many different calculations to be carried out in 
parallel. 



Computational circuits 

Although several different theoretical models are useful 
for abstract descriptions of computers, one of the most 
convenient approaches for describing how to build small 
computers is the circuit model, in which bits interact 
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FIG. 2; (a) A circuit to compute the exclusive-OR (xor) func- 
tion, z — X XOR y, using only nand and NOT gates. This is 
not the best such circuit; simpler circuits are known, but this 
one preserves the basic structure seen in figure hi. (b) A circuit 
to implement a NOT gate, y — NOT x using a nand gate. By 
combining circuits (a) and (b) it is possible to implement xOR 
using only nand gates. Any other function may be computed 
in a similar fashion, and so nand is a universal gate. Note, 
however, that both circuits use implicit CLONE gates, while 
(a) also uses one implicit SWAP gate. 



through gates which implement Boolean logic operations. 
One traditional set of classical gates is the one bit not 
gate together with two different two bit gates, and and 
OR. These three gates are said to form an adequate set, 
in that any desired logic operation can be performed by 
building an appropriate circuit (see figure m for an exam- 
ple) using some combination of these three gates. 

In fact it is not necessary to use all these gates: they 
can themselves all be obtained using combinations of 
NAND gates (figure 0), and so the nand gate is univer- 
sal for classical computation. It is, however, necessary 



to proceed with some caution, as several other "imphcit" 
gates are also required, such as the clone gate, which 
makes a copy of a bit, and the swap gate, which allows 
two wires to cross one another. 

This description works well with conventional comput- 
ers, in which bit states are represented by voltages ap- 
plied to wires, but it cannot be used with atomic com- 
puters. Atomic computers represent bit values using the 
quantum states of atomic systems, so a logic gate can nei- 
ther create nor destroy bits; thus logic gates such as and, 
which have two input bits and only one output bit, are 
immediately ruled out. Similarly, atomic systems evolve 
under a series of unitary transformations, which corre- 
spond to reversible operations, while many of the opera- 
tions described above are clearly not reversible. In order 
to build atomic computers, therefore, it is necessary to 
use a different approach, using only reversible logic op- 
erations. 



y — (B — y' 



FIG. 3: The controlled-NOT gate, which plays a central role 
in reversible and quantum computation. The target bit is 
marked by a ffi symbol, indicating the close relationship with 
binary arithmetic modulo two; the control bit is marked by a 
dot and a vertical control line (the dot is important in draw- 
ings of larger circuits where a control line may have to cross 
lines representing other bits). 
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Reversible computation 

The theory of reversible computation [|l5|, [2^, ^ has 
been studied extensively and, perhaps surprisingly, it is 
simple to perform any logic operation in a reversible man- 
ner. The only irreversible operation required when per- 
forming a computation |2£ 



30 is the preparation of a 



well defined initial state, usually taken as having all bits 
in state 0; after this initialisation process the computa- 
tion can be performed entirely reversibly. 

The basic approach needed to achieve reversible logic 
can be summarised in two simple rules. First, any logical 
inputs to a gate must be preserved in the outputs; this 
is most simply achieved by copying them to output bits 
without change. Secondly, the output of the gate (here 
assumed for simplicity to comprise a single bit) must be 
combined in a reversible fashion with an additional aux- 
iliary bit, for example by adding the two bits together 
using binary arithmetic modulo two. 

Binary arithmetic modulo two, usually indicated by 
the symbol ©, has the useful properties that x ® = 
Q<S)x = X and a;©l = l©a; = NOT(a;), while x®x = Q. 
A simple example of a reversible logic gate based on this 
approach is the controlled-NOT gate (see figure ||), which 
has two inputs, x and y and two corresponding outputs 
x' and y' . The first input bit is copied to its output bit, 
so x' — X, and is also combined with the second input 
bit to give y' — x ® y. Thus a not gate is applied to 
the second bit (the target bit) if and only if the first bit 
(the control bit) is in state 1. Note that a controlled-NOT 
gate is completely reversible; indeed it can be reversed by 
simply applying the same gate again (that is, controlled- 
NOT is its own inverse). Furthermore, as x(By = x XOR y 
the controlled-NOT gate is just a reversible XOR gate. 

A more complex example is provided by the controUed- 
controUed-NOT gate, often called the Toffoli gate (figure 



FIG. 4: The Toffoli gate, which gives x' 

z' — z® {x AND y). 



X, y 



y, and 



0), which plays a central role in reversible logic. This gate 
has three inputs, x, y and z, and three corresponding 
outputs, x' , y' and z'. The first two input bits, which are 
the logical inputs, are copied unchanged, so that x' = x 
and y' ^ y. A not gate is then applied to the third bit 
if and only if both x and y are in state 1; hence z' = 
z (B {x AND y). Thus this gate can be considered as a 
reversible equivalent to the conventional AND and nand 
gates: to evaluate x and y just use a Toffoli gate with 
z = 0, while for a nand gate set z — 1. 

The combination of the Toffoli gate with the 
controlled-NOT and simple not gates forms an adequate 
set, in that it is possible to build any other gate using only 
a combination of these gates (in particular controlled- 
NOT gates can also be used to build the two implicit 
gates, CLONE and SWAP, as shown in figure ||). Indeed, 
the Toffoli gate is universal in its own right, as a Toffoli 
gate can be easily converted into a controlled-NOT gate 
by setting a; = 1, and to a simple not gate by setting 
x = y = 1. 
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FIG. 5: (a) The controlled-NOT gate can be used to reversibly 
clone a bit; (b) three controlled-NOT gates implement a SWAP 
gate. 
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FIG. 6: The /-controUed-NOT gate, for reversible function 
evaluation: a' — a and b' = b (B f{a)- 



Reversible function evaluation 

Central to reversible computation is the idea of re- 
versible function evaluation. I will initially assume that 
the function has a single bit as both input and out- 
put, but the generalisation to more complex functions 
is straightforward. This can be achieved by constructing 
a circuit with two inputs and two outputs, as shown in 
figure H (an /-controUed-NOx) and setting b — 0. The 
two values of the function, /(O) and /(I), can then be 
evaluated by setting a — and a = 1 respectively. The 
/-controUed-NOT gate can itself be built out of simpler 
gates, such as those described above; in most cases this 
will also require a number of ancilla bits to hold interme- 
diate results. For simplicity these ancilla bits are usually 
omitted from diagrams such as figure pi 



QUANTUM COMPUTATION 



states contributing to the superposition at the same time, 
so a qubit can simultaneously occupy two different logical 
states. 

Computational circuits are implemented within quan- 
tum computers by performing physical manipulations so 
that the computer evolves under a propagator which im- 
plements the desired unitary transformation. Just as cir- 
cuits can be built up from gates these propagators can 
be assembled from simpler elements, and so these propa- 
gators are often referred to as circuits, even though their 
physical implementation may bear little resemblance to 
conventional electrical circuits. As before, this abstract 
model allows quantum computers to be described in a 
device-independent fashion. 

As discussed below, it is possible to construct any 
quantum circuit by combining a small number of sim- 
ple propagators, usually called gates. These gates only 
affect one or two qubits at a time, and so it is per- 
fectly practical to describe their propagators explicitly, 
for example as a matrix. A matrix description clearly 
depends on the choice of basis set, but the usual choice 
is to work in the computational basis, in which the ba- 
sis vectors correspond with the different logical states of 
the computer; thus for a two qubit gate the basis set is 
{|00), |01), 1 10), 1 11)}. In many proposed physical imple- 
mentations of quantum computers, such as NMR, this 
basis set is also the natural basis for the system, as the 
basis states are eigenstates of the background Hamilto- 
nian. 



We now have all the basic elements needed to describe 
how quantum computers could be used to extend our 
computational abilities. While large scale quantum al- 
gorithms, such as Shor's algorithm, are too complex to 
describe here, the basic ideas are relatively simple. 

As described above, a quantum mechanical two-level 
system can be used to build a reversible classical com- 
puter by using the two eigenstates of the system to rep- 
resent bits in logical states and 1. For example, the two 
Zeeman levels of a spin-^ nucleus in a magnetic field, \a) 
and 1/3), would be suitable for this purpose. For simplicity 
the two states are usually denoted |0) and |1), allowing 
quantum computation to be described without reference 
to any particular implementation; this choice of basis set 
is called the computational basis. The system will not be 
confined to these two eigenstates, however, but can also 
be found in superpositions, such as 



|0) + |i) 

V2 



(1) 



(the \/2 term is necessary to ensure that the state is 
normalised). For this reason a quantum mechanical two 
level system has much more freedom than a classical bit, 
and so is called a quantum bit, or qubit for short. A qubit 
in a superposition state is (in some sense) in both of the 



Quantum parallelism 

The central feature underlying all quantum algorithms 
is the idea of quantum parallelism, which in turns stems 
from the ability of quantum systems to be found in super- 
position states. Consider once again the reversible circuit 
for function evaluation (figure y). This can be achieved 
by constructing a propagator, Uf, applied to two qubits, 
which performs 

\a)\b)^\a)\b®fia)) (2) 

and setting \b) — |0). The two values of the function, 
/(O) and /(I), can then be evaluated by setting \a) = |0) 
and |a) = |1) as before. Now consider the effect of apply- 
ing this circuit when the first qubit is in a superposition 
described by equation |l| and the second qubit is in state 
|0) . This can be easily calculated as quantum mechanics 
is linear, and so the effect of applying a gate to a super- 
position is a superposition of the results of applying the 
gate to the two eigenstates. Hence the result is 

M±M|o) ^ 10)1/(0)) + |l)l/(l)) .3^ 

^/2 V2 

Thus the quantum computer has simultaneously evalu- 
ated the values of /(O) and /(I). 



When applied to more complex systems, quantum par- 
allelism has potentially enormous power. Consider a 
function for which the input is described by n bits, so 
that there are 2" possible inputs (since n bits can be 
used to describe any integer between and 2" — 1); a 
quantum computer using n qubits as inputs can evaluate 
the function over all of these 2" inputs in one step. In ef- 
fect a quantum computer with n input qubits appears to 
have the computational power of 2" classical computers 
acting in parallel. Unfortunately it is not always pos- 
sible to use this quantum parallelism in any useful way. 
Performing parallel function evaluation over n qubits will 
result in a state of the form 



2»/2 



(4) 



which is a superposition of the 2" possible inputs, each 
entangled with its own function value. If any attempt 
is made to measure the state of the system, then the 
superposition will collapse into one of its component val- 
ues, |r)|/(r)), where the value of r is chosen at random. 
Thus even though it seems possible to evaluate the func- 
tion over its 2" input values, it is only possible to obtain 
one of these values. 

It is, however, sometimes possible to obtain useful in- 
formation in a more subtle way. In some cases, the an- 
swer of interest does not depend on specific values, f{r), 
but only on global properties of the function itself. This 
is the basis of both Deutsch's toy algorithm and Shor's 
quantum factoring algorithm. 



Deutsch's algorithm 

Deutsch's algorithm [g3[ |^ [2^ was the first quantum 
algorithm to be discovered, and is one of the few quantum 
algorithms simple enough to describe here. The problem 
can be described in terms of function evaluation, but a 
more concrete picture can be obtained by thinking about 
coins. Normal coins have two different faces, convention- 
ally called heads and tails, but fake coins can be obtained 
which have the same pattern on both faces. 

Consider an unknown coin, which could be either a 
normal coin or a fake coin. In order to determine which 
type it is, it would seem necessary to look at both sides to 
find out whether they showed heads or tails, and then see 
whether these two results were the same (a false coin) or 
different (a true coin). With a quantum device, however, 
it would be possible to look at both sides simultaneously, 
and thus determine whether the coin was normal or fake 
in a single glance. 

The trick lies in abandoning any attempt to determine 
the pattern shown on either side of the coin; instead one 
must simply ask whether the two faces are the same or 
different. This is a property not of the individual faces 



of the coin, but of the whole coin, and thus may be ex- 
tracted from a state of the kind described by equation ||. 
A more detailed explanation of this approach is given in 
Section . 



Quantum gates 

Just as classical reversible computations can be per- 
formed using circuits built up out of reversible gates, 
quantum circuits can be constructed using quantum 
Unlike classical circuits, however, quantum 
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gates 

circuits can include gates which generate and analyse 

qubits which are in superpositions of states. 

One such gate is the single qubit Hadamard gate, H, 
which implements the transformations 



|0) ^ (|0) + |1))/a/2 



ID 



(|0)-|l))/y2 



(5) 



(6) 



As discussed below, this is closely related to a 90° 
pulse, but differs in some subtle ways; in particular the 
Hadamard is its own inverse. 

The Hadamard gate is useful as it takes a qubit in an 
eigenstate to a uniform superposition of states. By anal- 
ogy one can define multi-qubit Hadamard gates which 
take a quantum register into a uniform superposition of 
all its possible values: 



100) ^ (100) + 101) + |10) + |ll))/2. 



(7) 



This can be easily achieved by applying a one qubit 
Hadamard to each qubit. 

Another important property of the Hadamard gate can 
be seen by examining the right hand sides of equations |5| 
and y. These differ only by the presence of a minus sign, 
so the two superpositions differ only by the phase with 
which the two eigenstates are combined. The ability of 
the Hadamard gate to convert such phase differences into 
different eigenstates plays a central role in many quantum 
algorithms. 



Universahty of quantum gates 

The gates we have seen so far are enough to explain 
some simple quantum algorithms; for example, Deutsch's 
algorithm can be described using only Hadamard gates 
and reversible function evaluation gates. It is, however, 
useful to consider other more general quantum gates; in- 
deed as any unitary transformation can be considered 
as a quantum gate, we may need to consider an infinite 
number of gates. 

Within classical models of computation (both re- 
versible and irreversible) it is possible to construct any 
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desired gate by combining copies of a small number of 
simple gates. A similar situation applies in quantum 
computation, but in this case there are an infinite number 
of possible gates (even if we restrict ourselves to single 
qubit gates, any rotation around any axis constitutes a 
valid one qubit gate). Clearly it cannot be possible to 
construct an infinite number of different gates by com- 
bining a finite number of simpler gates, but it is possible 
to simulate any gate to any desired accuracy |^, |3J, 
which is good enough. Perhaps surprisingly there exists 
a very large number of two qubit gates which are uni- 
versal in this restricted sense, in that it is possible to 
simulate any desired gate (that is, any unitary transfor- 
mation) using only one of these universal gates together 
with its twin, obtained by swapping the roles of the two 
qubits. Indeed, mathematically speaking, almost all two 
qubit gates are universal |3J, |3^, 36 . 

While mathematically interesting, this result is of little 
immediate practical implication for most possible imple- 
mentations of a quantum computer, as it is usually more 
sensible to use a larger and more convenient set of gates. 
As one qubit gates are usually much simpler to perform 
than gates involving two or more qubits, it is often rea- 
sonable to assume that any one qubit gate (or, at least 
a reasonable approximation to it) is available. The com- 
bination of this set of one qubit gates with any single 
non-trivial two qubit gate, such as the controUed-NOT 
gate forms an adequate set [ p3[ , from which any other 
gate may be built with relative ease. 



BUILDING NMR QUANTUM COMPUTERS 

While it would in principle be possible to use a wide 
range of different approaches to build a quantum com- 
puter, all the main proposals to date 0, 0] have used 
broadly similar approaches, based on the quantum cir- 
cuit model outlined above. This model contains five ma- 
jor components, each of which must be implemented in 
order to construct a working computer [p7| . Four central 
components can all be implemented within NMR systems 
as described below, while the fifth component, error cor- 
rection, is discussed in Section . 



Qubits 

The first of these requirements, a set of qubits, ap- 
pears easy to achieve, as the two spin states of spin- 2 
nuclei in a magnetic field provide a natural implemen- 
tation. However, one important feature which distin- 
guishes NMR quantum computers from other suggested 
implementations is that NMR studies not a single iso- 
lated quantum system, but rather a very large number 
(effectively an ensemble) of such systems. Thus an NMR 
quantum computer is actually an ensemble of indistin- 



guishable computers, one on each molecule in the NMR 
sample. This has a number of subtle and important con- 
sequences as discussed below. 



Logic gates 

In order to perform an arbitrary computation it is 
necessary to implement arbitrary quantum logic circuits. 
This can be achieved as long as it is possible to imple- 
ment an adequate set of gates, which can be combined 
together to implement any other desired gate. While 
many different sets of gates are possible, a simple ap- 
proach is to implement the set of all possible one qubit 
gates, together with one or more non-trivial two qubit 
gates ||. 

One qubit gates correspond to rotations of a spin 
within its own one-spin Hilbert space, which can be read- 
ily achieved using RF fields. Note that it is necessary to 
apply these rotations selectively to individual qubits. In 
most other suggested implementations of quantum com- 
putation H, H this is easily achieved using some type of 
spatial localisation: the physical objects implementing 
the qubits are located at well defined and distinct loca- 
tions in space. This approach is not possible in NMR, as 
each qubit is implemented using an ensemble of nuclei, 
each of which is located at a different place in the NMR 
sample, and all of which are undergoing rapid motion. 
Instead different qubits are implemented using different 
nuclei in the same molecule, and they are distinguished 
using the different resonance frequencies of each nucleus. 

Two qubit gates, such as the controUed-NOT gate, are 
more complicated as they involve conditional evolution 
(that is, the evolution of one spin must depend on the 
state of the other spin), and thus require some interac- 
tion between the two qubits. The J-coupling in NMR 
is well suited to this purpose. Note that all the differ- 
ent nuclei making up an NMR quantum computer must 
participate in a single coupling network. It is not nec- 
essary (or even desirable) that all the nuclei are directly 
coupled together, but they must be connected, directly 
or indirectly, by some chain of resolved couplings. Since 
J-coupling only occurs within a molecule, and does not 
connect different molecules, we can treat an ensemble of 
molecules as an ensemble of identical mutually isolated 
computers. 



Initialisation 

Quantum logic gates transform qubits from one state 
to another, but this is only useful if the qubits start off in 
some well defined initial state. In practice it is sufficient 
to have some method for reaching any one initial state, 
and the obvious choice is to have all the qubits in the |0) 



state, corresponding to a CLEAR operation. Any other 
desired starting state can then be easily obtained. 

When, as for NMR, the computational basis coincides 
with the natural basis of the quantum system it should 
in principle be easy to implement CLEAR as it takes the 
quantum computer to its energetic ground state, and this 
can be achieved by some cooling process. Unfortunately 
this approach is not practical in NMR as the Zeeman en- 
ergy gap is small compared with the Boltzman energy at 
any reasonable temperature; thus at room temperature 
the population of all the states will be almost equal, with 
only small deviations (around one part in 10'*) from the 
average. Techniques for enhancing spin polarization [p8|, 
such as optical pumping [^ ^, ^ , and the use of para- 
hydrogen ^, ^ Q allow this deviation to be increased, 
but with the exception of optically pumped noble gases 
it has so far proved impossible to even approach a pure 
ground state system. 

This apparent inability to implement the clear op- 
eration led to NMR being rejected as a practical tech- 
nology for implementing quantum computers. Recently, 
however, it was realised |q| that this conclusion was over 
hasty, as with an ensemble quantum computer it is not 
actually necessary to produce a pure ground state; in- 
stead it suffices to produce a state which behaves in the 
same manner as the pure ground state. This point can 
be clarified by considering the density matrix describing 
a single isolated spin-half nucleus. This exhibits nearly 
equal populations for the two eigenstates, but with a 
slight excess in the (low energy) |0) state compared with 
the (slightly higher energy) |1) state. No NMR signal will 
be observed from the equal populations, as the signals 
from different molecules will cancel out, but a small signal 
can be seen which arises from the deviations away from 
the average. Thus, ignoring questions of signal intensity, 
for a single isolated nucleus the thermodynamic equilib- 
rium state is indistinguishable from a pure |0) state. 

States of this kind are often called pseudo-pure states, 
or effective pure states [g, |[ |o| . Unfortunately the sim- 
ple approach outlined above does not work for larger spin 
systems, as the pattern of population deviations is more 
complicated, and does not have the desired form. Sev- 
eral different techniques have, however, been developed 
to tackle this problem. 



Readout 

The last stage in any quantum computation is to char- 
acterise the final state of the system, so that the result of 
the computation may be read out. Just as for initialisa- 
tion, a range of different approaches have been used, but 
all these approaches combine two major elements. For 
simplicity I will assume that the computation ends with 
the result qubits in eigenstates; thus it is only necessary 
to determine whether a given qubit is in (the pseudo- 



pure) state |0) or |1). 

The simplest approach is to apply a 90° pulse to the 
corresponding spin, and observe the NMR spectrum Q . 
Since |0) corresponds to the ground state, a qubit in |0) 
will give rise to an absorption line; correspondingly a 
qubit in state |1) will give an emissive signal. It is, of 
course, necessary to acquire some sort of reference signal, 
in order to distinguish between these two extremes, but 
this can be easily achieved by acquiring the spectrum of 
the pseudo-pure initial state. 

The second major approach p2| is to determine the 
state of one qubit by analysing the multiplet structure 
within the spectrum of a neighbouring spin. If several 
spins are coupled together, then individual lines within 
a multiplet can be assigned to specific states of these 
neighbours. Thus, the spectrum of one spin can give 
information on the states of several different qubits. 



Some two spin systems 

While a number of different systems have been used 
to build small NMR quantum computers, all their ma- 
jor features can be explored using two different two-qubit 
systems which were used in the earliest demonstrations 
of NMR quantum computation [^l| |l^. The most im- 
portant difference between these systems is that one uses 
a homonuclear two-spin system, while the other is het- 
eronuclear. 

The first example system uses the two ^H nuclei of 
partially deuterated cytosine in D2O (see figure 0). As 
this system is homonuclear it is possible to excite both 
nuclei with a single hard pulse, and to observe both nu- 
clei in the same spectrum. Another more subtle advan- 
tage is that the pattern of Boltzmann populations is sim- 
pler in homonuclear systems than in their heteronuclear 
counterparts. There are, however, two significant disad- 
vantages of such as system. Firstly the two *H multi- 
plets have relatively similar frequencies, as they lie only 
about 1.51 ppm apart, and thus it is necessary to use 
soft frequency selective pulses |^ (or sequences of hard 
pulses and delays with equivalent effects) in order to ad- 
dress the spins individually. Secondly, the J-coupling 
between the two spins is relatively small (about 7 Hz), 
and so controlled gates take a fairly long time to im- 
plement. It would, of course, be possible to choose a 
different molecule, in which the chemical shift difference 
or J-coupling was larger, but it is difficult to improve 
one without making the other worse. While it is unlikely 
that cytosine is the absolutely optimal choice, no other 
homonuclear *H system would be very much better. 

The heteronuclear alternative is probably the most 
widely used two qubit NMR system. It is based on the 
^H and *'^C nuclei in ^^C-labeled chloroform. This has 
the huge advantage that it is possible to separately ex- 
cite the two spins using hard pulses, rendering selective 
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FIG. 7: The structure of partially deuterated cytosine ob- 
tained by dissolving cytosine in D2O; the three protons bound 
to nitrogen nuclei exchange with solvent deuterons, leaving 
two ^H nuclei as an isolated two spin system (all other nuclei 
can be ignored). 



excitation essentially trivial. Furthermore, the relatively 
large size of the J-coupling allows two qubit gates to be 
performed much more rapidly than in homonuclear sys- 
tems. In this heteronuclear system it is not possible to 
acquire signals from both spins simultaneously, but this 
is not a major problem as it is possible to determine the 
states of both spins by examining either the ^H or the 
^^C spectrum. Similarly, the complex pattern of popula- 
tions over the four energy levels of this system does not 
fit with the original scheme for generating pseudo-pure 
states; however, some more modern schemes are in fact 
simpler to implement in heteronuclear systems. 

Considering all these issues together, it is not easy to 
say whether it is better to use homonuclear or heteronu- 
clear systems to implement two qubit NMR quantum 
computers: heteronuclear systems are perhaps simpler to 
work with, but homonuclear systems give more elegant 
results. With larger spin systems the issues become even 
more complex, and a wide range of options have been ex- 
plored. It is clear, however, that the simplest approach 
of using a fully heteronuclear spin system is unlikely to 
be practical beyond five qubit systems, as there are only 
5 "obvious" spin-half nuclei which can be used (^H, ^'^C, 
^^N, ^^F and '^^P). In practice NMR quantum computers 
with more than three qubits are likely to include two or 
more spins of the same nuclear species; it is, therefore, 
essential to consider how computation can be performed 
in homonuclear systems. 



Scaling the system up 

The requirements outlined above are adequate for 
building small quantum computers, suitable for simple 
demonstrations of quantum information processing. If, 
however, one wishes to build a large scale quantum com- 
puter, suitable for performing interesting computations, 
then it is necessary to consider whether the approaches 
used are limited to such small systems, or whether (and 
if so, how) they can be scaled up. A fifth requirement for 
practical quantum computation B^ , the implementation 



of fault-tolerant quantum error correction, is described 
in Section . 

This is an important practical question, but not one 
which will be addressed in detail here. The problems of 
scaling up NMR quantum computers are formidable, and 
have been well described elsewhere [38, 4^, ^. Most au- 
thors now agree that NMR approaches are likely to be 
limited to computers containing 10-20 qubits; this is sig- 
nificantly smaller than estimates of the size required to 
perform useful computations (50-300 qubits). Further- 
more the apparent inability of NMR systems to perform 
efficient quantum error correction rules out their use for 
many types of problem. 

The fundamental difficulties involved in scaling up cur- 
rent NMR quantum computers to large sizes have led 
some authors to suggest that this approach does not ac- 
tually implement real quantum computation at all. This 
is a quite subtle question which will be discussed further 
in Section below. 



QUBITS AND NMR SPIN STATES 

Traditional designs for quantum computers comprise 
a number of two-level systems which interact with one 
another and have some specific interaction with the out- 
side world, through which they can be monitored and 
controlled, but are otherwise isolated. NMR systems are 
rather different: a typical NMR sample comprises not 
one spin-system, but a very large number of copies, one 
from each molecule in the sample, effectively forming an 
ensemble of copies. Traditional quantum computers are 
usually described using Dirac's bra(c)ket notation [[l7[ , 
but NMR systems are better described using density ma- 
trices, usually written in the product operator basis pq , 
which has a number of important consequences. It is pos- 
sible to draw close analogies between the states of tradi- 
tional quantum computers and those used in descriptions 
of NMR systems ^8| , but it is necessary to proceed with 
caution. 



One qubit states 

A single qubit can be in either of its two eigenstates, 
|0) and |1), or in some linear superposition of them. Such 
a state is most conveniently written as a column vector 
in Hilbert space, for example 



(8) 



NMR quantum computers cannot be properly described 
in this way; instead they must be described using the 
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corresponding density matrix 

p^\4'){i>\-- 



CqCo CjCo 

CqCi clci 



(9) 



which can then be decomposed as a sum of the four Pauh 
basis states or their product operator equivalents, ^E, I^, 
ly, and J^. 

Consider first the eigenstates, |0) and |1), which corre- 
spond to the density matrices 



and 



|1)(1| 



1 





1 



k^E- 



\E- 



(10) 



(11) 



respectively. As all NMR observables are traceless, mul- 
tiples of the unit matrix can be added to density matrices 
at will, and so as far as any NMR experiment is concerned 
the density matrix I^ is equivalent to |0), while — /^ is 
equivalent to |1). In the language introduced above, Iz 
and —Iz are pseudo-pure states, corresponding to |0) and 
|1) respectively. This approach cannot, however, be ex- 
tended to larger spin systems without modifications. 

Next consider superpositions, such as (|0) -I- |l))/\/2, 
with its corresponding density matrix 




\E 



/x 



(12) 



As before multiples of the unit matrix can be ignored, and 
so (|0)-|-|1))/a/2 is equivalent to Ix- Similarly |0)-|-i|l) is 
equivalent to ly, while |0) — |1) is equivalent to —Ix- Just 
as single qubit eigenstates are closely related to one spin 
magnetizations, their superpositions are closely related 
to one spin coherences. 



Two qubit states 

While there is a simple relationship between qubit 
states and NMR states for a single qubit (a one spin sys- 
tem), this relationship is more complicated in systems 
with two or more qubits |}48| . Typically quantum al- 
gorithms start with all qubits in state |0), which for a 
two-qubit computer is the state |00). The corresponding 
density matrix 



|00)(00| = 



/I 0\ 





\0 0/ 



(13) 



is not the same as the thermal equilibrium density matrix 

^10 Q\ 



yo -ij 



Iz+Sz^ 



(14) 



The ideal density matrix (Eq. |l^) can, however, be de- 
composed as the sum of four product operators: 



1 00) (00 1 = 



hE 



Iz+Sz+ 2IzSz 



(15) 



and this sum (ignoring multiples of the unit matrix as 
usual) can be assembled using conventional NMR tech- 
niques, as described below. 

Superpositions can be treated in nruch the same way, 
but they are not directly related to NMR coherences in 
any very simple way. For example consider the state 
(|00) + |01))/\/2, in which the first spin is in state |0), 
while the second spin is in a superposition of states. 
The corresponding density matrix can be decomposed 
directly: 



/^ 



0\ 





Vo oy 
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1 1 

2 2 



^{^E + Iz + Sx + 2IzSx), (16) 



but there is a more subtle approach. Note that (|00) + 
|01))/V2 can be written as a product of single qubit 
states 



|00) + |01) |0)(|0)-f|l)) 



V2 



V2 



(17) 



and so the corresponding density matrix can also be de- 



composed as a direct product of equations nO and 12 
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^Sx) 



(18) 



Unlike the single qubit case, a simple superposition does 
not correspond directly to an NMR coherence, but in- 
stead to a complex mixture of coherences and popula- 
tions. It is, however, rarely necessary to worry about 
this, as such states can be easily obtained from states 
like Eq. |l|. 

Finally consider superpositions of the form (|00) -|- 
|ll))/v2, which cannot be broken down into a product 
of one qubit states (such states are said to be entangled). 
As they cannot be factored it is necessary to decompose 
the corresponding density matrices directly. In this case 
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Zil^tJ^; ~l ^Ix^X 



2IySy) 



(19) 

which is a mixture of longitudinal two-spin order and 
DQx double quantum coherence. 
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NMR LOGIC GATES 

After the rather abstract discussions above, we now 
turn to the details of methods by which quantum logic 
gates can be (and have been) implemented within NMR. 



One qubit gates 

Many one qubit logic gates can be implemented di- 
rectly. For example, a simple not gate, which intercon- 
verts |0) and |1), can be implemented as a 180° rota- 
tion Q. Rotations about axes in the xy- plane can be 
achieved using RF pulses, while rotations about the z- 
axis can be accomplished either by using periods of free 
precession under the Zeeman Hamiltonian, or by compos- 
ite z-pulses Q. This does not, however, cover the full 
range of gates which may be desired, as some of these 
correspond to rotations about tilted axes. 

An obvious (and important) example is the Hadamard 
gate. While this superficially resembles a 90° pulse, this 
resemblance is misleading, as the Hadamard gate is its 
own inverse. Clearly the Hadamard must correspond to 
a 180° rotation, and a little thought reveals that this 
rotation occurs around an axis tilted at 45° within the 
xz-plane. This could be achieved directly by using off- 
resonance excitation, but this has a number of practical 
difficulties. Alternatively it can be implemented using 
a composite pulse sequence, such as 45°-180°-45f_y; as 
180°-45j_y can be replaced by 45^^-180° this three pulse 
sequence may be simplified to the two pulse sequence 
90°-180° or 180°-90^j^. 

In fact, when implementing quantum algorithms on 
NMR quantum computers it is rarely necessary or desir- 
able to use a Hadamard gate, as it can generally be re- 
placed by the NMR pseudo-Hadamard gate, a 90° pulse 
p8[ . As this gate is not self- inverse it is usually neces- 
sary to replace pairs of Hadamard gates by one pseudo- 
Hadamard and one inverse pseudo Hadamard (901y) 
gate. This is a simple example of a general rule in experi- 
mental implementations of quantum computation: rather 
than directly implementing the gates commonly used in 
theoretical descriptions, it is better to use simpler gates 
which are broadly functionally equivalent to them. 



(a) 



Controlled-NOT gates 

This approach is also applicable to the implementation 
of controlled two-qubit gates. While it is perfectly possi- 
ble to implement a controUed-NOT gate, this is not nec- 
essarily the most sensible approach. The controUed-NOT 
gate can itself be assembled from simpler basic gates, and 
it may be more sensible to use these basic gates directly. 

A natural way to implement a controUed-NOT gate is 
to use a three gate circuit, as shown in figure |3(a). The 
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FIG. 8: (a) A circuit for implementing a controlIed-NOT gate; 
the two outer gates are one qubit Hadamard gates, while the 
central gate is a controlled vr-phase-shift gate (see text for de- 
tails), (b) A circuit for implementing a controIled-NOT gate on 
an NMR quantum computer; the one qubit Hadamard gates 
are replaced by pseudo-Hadamard (h) and inverse pseudo- 
Hadamard (h~^) gates. 



two boxes marked H are one qubit Hadamard gates, and 
the central gate (two circles connected by a control line) 
is a two qubit controlled vr phase-shift gate. This gate 
performs the transformation 



11)11) 



-11)11) 



(20) 



while leaving all other states unchanged, and so is de- 
scribed by the matrix 



/I \ 
10 
10 

\0 -1/ 



(21) 



Unlike the controUed-NOT gate this phase shift gate 
is symmetric; it is not meaningful to ask which qubit 
the phase shift was applied to. Note that the other 
controUed-NOT gate, in which the roles of control and 
target qubit are reversed, can be constructed by simply 
moving the two Hadamard gates to the upper line. 

When implementing this circuit on an NMR quantum 
computer it is preferable to avoid using Hadamard gates, 
as these are difficult to implement. Instead these two 
gates are replaced by an inverse pseudo-Hadamard (a 
90_y pulse) and a pseudo-Hadamard gate (a 90j^ pulse) 
respectively, as shown in figure |8|(b) . 

To see how to implement the controlled phase shift 
gate it is best to break down the propagator, equation 
|2l| , using the product operator basis set. As this propa- 
gator is diagonal, it must arise from evolution under the 
diagonal operators I^, 5^, 21 ^Sz and (for completeness) 
kE. It can be decomposed in two ways. 



TT = exp [-i X 7r/2 {hE -1^-3^ + 21^8^ 



(22) 



TT = exp [-i X 7r/2 {-^E + 1^ + 3^- 21^3^)] ■ (23) 

The choice between these two decompositions is simply 
a matter of experimental convenience. 

The use of 5 -E in these equations might appear to give 
rise to difficulties, as this is not normally considered as a 
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product operator which can arise in NMR Hamiltonians. 
In fact it is of no significance at all, as the effect of the iE 
term is simply to impose a global phase shift. Such global 
phase shifts have no physical meaning, and cannot be de- 
tected; note that these global phase shifts have no effect 
on density matrix (or product operator) descriptions of 
the spin system. Physically this corresponds to the fact 
that there is no absolute zero against which energies may 
be measured. 

Implementing these propagators (equations 23 and p3) 
is fairly straightforward. The pure spin-spin coupling 
term (21 ^Sz) can be generated using conventional spin- 
echo techniques, while the two Zeeman Hamiltonians {Iz 
and Sz) can be achieved by appropriately timed periods 
of free precession, by the use of composite z-pulses, or 
most simply by just rotating the RF reference frame. As 
the three terms all commute, they need not be applied 
simultaneously, but can be applied in any order. When 
the individual elements making up the propagator are 
combined, it is frequently possible to combine or cancel 
individual pulses, thus simplifying the whole sequence. 
This results in a wide variety of possible pulse sequences, 
and the choice among them is largely a matter of taste. 
For example, one possible sequence for implementing w 



1 
4J 



ISO^ 



1 
4J 



90^ - 90_ 



90. 



(24) 



where all pulses are applied to both spins. 

While the possible pulse sequences differ in detail they 
have one feature in common: an evolution time of 1/2 J 
(occasionally 3/2 J) during which the spins evolve under 
the spin-spin coupling, so that the antiphase condition is 
achieved. This is, of course, the central feature of coher- 
ence transfer sequences, such as INEPT, indicating the 
close relationship between controlled two qubit gates and 
coherence transfer. 



Other two qubit gates 



with (f> — n/n. Clearly (p can be constructed in much the 
same way as tt, equation k3. 



Gates in larger spin systems 

The approaches described above can be easily imple- 
mented in two spin systems, allowing quantum computers 
with two qubits to be easily constructed. With larger 
spin systems, however, the process can become much 
more complicated [^. It is not possible simply to use 
pulse sequences designed for two spin systems, as it is 
necessary to consider the evolution of all the additional 
spins in the system. In particular it may be necessary to 
refocus the evolution of these spins under their chemical 
shift and J-coupling interactions. The simplest method 
is to nest spin echoes within one another, so that all the 
undesirable interactions are removed, but this naive ap- 
proach requires an exponentially large number of refocus- 
ing pulses (that is, the complexity of the pulse sequence 
doubles with every additional spin) . This problem can be 
overcome by using efficient refocusing sequences pO| , p^ , 
which allow refocusing to be achieved with quadratic 
overhead. 

It is, of course, rare to find a large spin system where 
all the couplings have significant size; in most cases long 
range couplings will be small enough to be neglected. 
This greatly simplifies the problem, both by reducing 
the number of couplings which have to be refocused, and 
by simplifying the echo sequences required |50, |2[. It 



might seem that it would be difficult to implement some 
logic gates in such a partially coupled spin system, as 
the necessary spin-spin couplings are missing. In fact 
this is not a problem, as long as every pair of spins is 
connected by some chain of couplings: quantum swap 
gates p3, M] can be used to move quantum information 
along this chain. 



While the controUcd-NOT gate is important it is not the 
only two-qubit gate worth considering: while it is possible 
to construct any desired gate using only controUed-NOT 
gates and one-qubit gates, it is usually more efficient to 
use a wider repertoire of basic gates. 

One simple and important example is the controlled 
square-root of not gate, which plays a central role in 
traditional constructions of the three-qubit Toffoli gate 
pl[ ; this is one member of a more general family of n 
roots of NOT. Such gates can be built in much the same 
way as controUed-NOT gates, except that the controlled 
TT phase-shift gate must be replaced by the more general 
transformation 
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(25) 



Multi qubit logic gates 

Multi qubit logic gates are gates, such as the Tof- 
foli gate, which perform controlled operations involving 
more than two qubits. Such gates can of course be im- 
plemented by constructing appropriate networks of one 
qubit and two-qubit gates pit] , but as the NMR Hamilto- 
nian can contain terms connecting multiple pairs of spins 
it should be possible to build some such gates directly, 
with a significant saving in pulse sequence complexity. 
This is indeed the case, and several interesting results 
have been obtained |5^, |5^ . 
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Single transition selective pulses 

An alternative approach for building controlled two 
qubit gates is to use ultra-soft selective pulses [Q, with 
excitation profiles so narrow that they pick out, for ex- 
ample, a single transition from within a doublet. This 
corresponds to only exciting the nucleus of interest when 
the neighbouring nucleus is in a certain state fST]] ; if the 
excitation corresponds to a 180° pulse then this provides 
a simple way of constructing a controUed-NOT gate [^ . 
The relationship between this approach and the (more 
common) multiple pulse sequence approach is analogous 
to that between the old fashioned selective population 
transfer experiment [ p9| and its more modern counter- 
part, INEPT @. 

One advantage of this approach (58| Q is that it is 
relatively easy to extend it to multi qubit gates such as 
the Toffoli gate. This can be achieved by using a selective 
pulse which affects one of the four transitions of a spin 
coupled to two neighbours. A corresponding disadvan- 
tage is that in this case constructing a simple controUed- 
NOT gate requires either a selective pulse which excites 
two of the four transitions in such a system or the applica- 
tion of two single transition selective pulses in sequence. 

The traditional Toffoli gate corresponds to inverting 
a qubit when two other qubits are in the state |1), but 
there is an entire family of related gates which effect an 
inversion for some given pattern of states. Each such gate 
corresponds to exciting a different transition in the mul- 
tiplet, and so the entire family of gates can be achieved 
directly. In practice, however, the central regions of a 
multiplet can become quite crowded, with many lines 
nearly overlapping, and it will be difficult to select a sin- 
gle transition. By contrast the two transitions at the 
extreme ends of the multiplet will always be relatively 
well separated from their nearest neighbours, and it is 
best to concentrate on these two frequencies. Other gates 
can then be constructed by surrounding these basic gates 
with NOT gates applied to the neighbouring spins, thus 
permuting the identities of the lines in the multiplet. 



Geometric phase-shift gates 

A third approach for implementing NMR quantum 
computation, based on the use of geometric phase-shift 
gates, has recently been described ^^ |6^. Like the 
conventional approach it relies on controlled phase-shift 
gates, but the phase shifts are generated using geometric 
phases p3, such as Berry's phase pa, rather than the 
more conventional dynamic phases. Berry phases have 
been demonstrated in a wide variety of systems pM| , in- 
cluding NM R ||66|, |67t and the closely related technique 
of NQR pi, |69|, |70(| , and can be used to implement con- 
trolled phase shift gates in NMR systems |6l, Bsl. This 



approach has few advantages for NMR quantum compu- 
tation, but may prove useful in other systems [|63[ . 



INITIALISATION AND NMR 

As it is impractical to cool down NMR spin systems to 
their ground state [|8[ ^, Q , initialisation of an NMR 
quantum computer in practice means assembling an ap- 
propriate pseudo-pure state. This approach is useful only 
if some practical procedure for assembling such states can 
be devised. 

For the simplest possible system (a single nucleus) the 
process is trivial, as the thermal equilibrium density ma- 
trix has the desired form, but with larger systems the 
situation is more complicated. The essential feature of a 
pseudo-pure state is that it has a diagonal density ma- 
trix in which the populations of all the spin states (the 
elements along the diagonal) are the same, with the ex- 
ception of one state (normally |0) = |000 ... 0)) which has 
a larger population. By contrast, at thermal equilibrium 
the spin state populations are distributed in accordance 
with the Boltzmann equation, and so exhibit a more com- 
plex variation. For a homonuclear two spin system the 
equilibrium density matrix (neglecting multiples of the 
identity matrix and an initial scaling factor) is Iz + Sz, 
while the desired pseudo-pure state is proportional to 
Iz + Sz + 2IzSz (equation ^. 

The original approach for assembling pseudo-pure 
states, developed by Cory et al. [g|, g|, uses conventional 
NMR techniques. Assembling such a mixture using pulse 
sequences and field gradients is a fairly straightforward, if 
somewhat unusual, NMR problem. This process is com- 
monly called spatial averaging, presumably a reference to 
the use of field gradients. 

A second early approach, suggested by Gershenfeld 
and Chuang Q, is to use a subset of the energy lev- 
els in a more complex spin system. For example, in a 
homonuclear three spin system it is possible to find a set 
of four energy levels which exhibit the pattern of popu- 
lations corresponding to the pseudo-pure state of a two 
spin system. This approach, often called logical labeling, 
is elegant in principle but complex to apply in practice, 
and has only rarely been experimentally demonstrated. 

More recently a variety of different approaches have 
been used, although these all combine elements of the 
two basic approaches above. The most popular tech- 
nique, usually called temporal averaging ]71| |, works by 
performing many different experiments, each with a dif- 
ferent initial state. For example, in a two qubit system, 
one might perform experiments starting from Iz, Sz, and 
2IzSz- If the spectra from these three experiments are 
then added together, the result is equivalent to a sin- 
gle experiment starting from a mixture of these states. 
Clearly temporal averaging and spatial averaging are re- 
lated in much the same way as coherence selection meth- 
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ods based on phase cycling and gradients. Finally, a new 
approach combines these methods in a cunning way, us- 
ing the analogy between multiple quantum coherence and 
so-called "cat" states to generate pseudo-pure states in a 
fairly efficient manner p2[. 



Spatial averaging 

The direct "spatial averaging" technique may be ex- 
emplified by the original sequence of Cory et al. M, m for 
constructing a pseudo-pure state in a two spin system: 
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2 0^ — 2 Oil 



45°/, 



-^Iz + ^Sz 



couple 



45°/-. 



\/2-^z V^-'-y 






V2-^z 



2 0z 
-^ZlxJz + 20z 



(26) 



11 z ^ 2-/a; + 2 2Ixbz + 2 0z + 2 2,Izbz 



crush I J lie I lore 

'*2-lz I 20z + 2^1ziJz 



where the sequence is described in product operator no- 
tation, "crush" indicates the application of a crush field 
gradient pulse, and "couple" indicates evolution under 
the scalar coupling for a time 1/2 J. Note that the two 
crush pulses must be applied along different axes, or with 
different strengths, to prevent undesired terms from be- 
ing refocused. 

Several alternative sequences for creating two spin 
pseudo-pure states have been developed; for example fz^ 



Iz+Sz 



45°(/:„+Sx) couple 30°(7_„+S'_a) 



{Iz +Sz+ 2IzSz 



(27) 



where zero quantum terms (which in a homonuclear spin 
system will survive the crush pulse) have been neglected. 
This scheme works well in heteronuclear spin systems, 
but in homonuclear systems it is necessary to use a more 
complex approach to deal with the zero quantum terms. 
These sequences can be generalised to larger spin sys- 
tems, but this process is quite complex. For this reason 
most work on larger spin systems has used temporal av- 
eraging techniques. Recently, however, Knill et al. Q 
have developed a general scheme based on cat states, 
which allows pulse sequences for any spin system to be 
developed. This approach is described below. 



Logical labeling 



Logical labeling |10 is most easily understood by ex- 



homonuclear three spin system: 

Iz + Sz+Rz^h {3, 1, 1^ -1^ 1, -1, -1, -3}, (28) 

where the braces indicate a diagonal matrix defined by 
listing its diagonal elements. While this matrix does not 
have the right form for a three spin pseudo-pure state it 
is possible to select out four levels (corresponding to the 
states |000), |011), |101) and |110)) which have the same 
population pattern as 



4 + 5, + 2/,5, = i {3, -1,-1,-1}, 



(29) 



and so this subset of levels can be used as a two spin 
pseudo-pure state. 

It would be possible to use these states directly, but 
this would greatly complicate subsequent logic opera- 
tions as there is no simple correspondence between these 
four states of the three spin system and the four basic 
states of a two spin system. Instead it is better to per- 
mute the populations of the various states, performing 
|001) ^^ |101) and |010) ^^ |110); these permutations can 
be achieved using controUed-NOT gates. At the end of 
this process the populations are given by 



i{3,-l, -1,-1, 1,1, 1,-3}, 



(30) 



amining the thermal equilibrium density matrix for a 



so that the states |000), |001), |010) and |011) are in a 
pseudo-pure state. Note that these four states all have 
the first spin in state |0), and so the first spin acts as an 
ancilla spin, labehng the "correct" subspace. 

Similar, but more complex, procedures can be used 
with larger spin systems Q. The overhead required 
is fairly small; that is the number of pseudo-pure spins 
which can be encoded in a spin system is only slightly 
smaller than the size of the system. However, while these 
results are elegant the complexity of implementing logical 
labeling means that experimental demonstrations have so 
far been confined to three spin systems ^^ [74[ . 



Temporal averaging 

As discussed previously, temporal averaging ||7l| bears 
much the same relationship to spatial averaging as phase 
cycling does to the use of gradients to select coherence 
transfer pathways. The name can, however, be used to 
cover a variety of different approaches. 

As described above (equation ^, a pseudo-pure state 
of a two spin system can be assembled as a mixture of 
three terms: Iz, Sz and 2IzSz- The simplest approach 
to temporal averaging is just to perform a computation 
starting from each of these states, and add the results 
together at the end. This is easily generalised to larger 
spin systems: for a system of n spins it is necessary to 
perform 2" — 1 separate experiments. In some simple ex- 
periments it is possible to show that only some of these 
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starting states will give an observable signal |Q , and so 
it is unnecessary to perform experiments starting in other 
states. This permits substantial experimental simplifica- 
tions, but it is not a general technique. 

A better approach is to use the original scheme of Knill 
et al. | |7l[ . The thermal equilibrium density matrix for a 
two spin system 



|0)- H 



4 + 5, = {1,0,0,-1} 



(31) 



(where the braces have the same meaning as before) can 
be easily converted into two other states, 



and 



/, +245, ={1,0,-1,0} 



5, + 245, = {1,-1,0,0}. 



(32) 



(33) 



These three states are related by simple permutations of 
the populations of the levels, which can be achieved using 
controUed-NOT gates. Adding together the three starting 
states gives 

2 (4 + 5, + 245,) = {3, -1, -1, -1} , (34) 

which is a pseudo-pure state. Adding together the spec- 
tra from computations started in these three states there- 
fore gives the spectrum which would be produced from a 
pseudo-pure state. 

Once again this process is easily generalised to larger 
spin systems. The most obvious approach is to average 
over the 2" — 1 cyclic permutations of the populations in 
an n spin system. This exhaustive averaging scheme is 
just as inefhcient as the naive approach outlined above, 
but Knill et al. [[7l| have shown that similar results can 
be achieved by averaging over much smaller numbers of 
states. 



The use of "cat" states 

The schemes described above are perfectly practical for 
small spin systems but are harder to use with larger sys- 
tems. Recently Knill et al. |Q have described a simple 
approach which works for spin systems of any size and 
which can be used with either the gradient (spatial aver- 
aging) or phase cycling (temporal averaging) approaches. 
Their method is based on the properties of "cat" states, 
named by analogy with Schrodingcr's Cat. An n qubit 
cat state is a superposition state of the form 



0± = (|OO...O)±|11...1))/V2, 



(35) 



so that either all the n qubits are in state |0), or all the 
qubits are in state |1). (In fact the relative phase of the 
two states contributing to the superposition can take any 
value between and 27r, but it is convenient to restrict 



-^ 



-^ 



FIG. 9: A network for creating the three qubit cat state (fi^ . 

ourselves to the two values and tt, giving rise to the fac- 
tor of ±1.) States of this form are said to be entangled, 
and play a central role in quantum information process- 
ing and experimental tests of quantum mechanics. The 
role of entanglement in NMR quantum computers will 
be explored in more detail below, but for the moment it 
is sufficient to note that cat states are closely related to 
(but not simply equivalent to) multiple quantum coher- 
ence Q. 

As discussed above (equation noh the two qubit cat 
state (j)^ = (|00) + |11))/V2 (commonly called a Bell 
state W^) is a mixture of DQ^ double quantum coher- 
ence and longitudinal two-spin order. Similarly the three 
qubit cat state 4>t (usually called a GHZ state |Q) is a 
mixture of ?)Qx triple quantum coherence and the three 
possible states of longitudinal two-spin order, and a gen- 
eral n qubit cat state will correspond to a mixture of n 
quantum coherence and ordered population states. Thus 
an n quantum filtration sequence is almost (though not 
quite) equivalent to selecting n qubit cat states. 

Cat states are easily prepared from pure states, using 
controUed-NOT gates. One possible network for a three 
qubit system is shown in figure |9|; networks for larger 
systems can be derived by analogy. Similarly by reversing 
this network cat states can be converted back into pure 
states. Thus, if it is possible to prepare an n qubit cat 
state, it should be possible to obtain a corresponding 
pure state. 

This suggests a simple scheme for preparing pseudo- 
pure states. If the network shown in figure is applied 
to a three spin system in its thermal equilibrium state, 
the resulting mixture will include a component of triple 
quantum coherence, and thus of the desired cat state. 
This component can be selected, either by phase cycling 
or by using gradient methods. Finally the network can be 
reversed to convert the cat state back into a pseudo-pure 
state. 

Unfortunately this does not quite have the desired ef- 
fect, as triple quantum coherence is not quite equiva- 
lent to the desired cat state; in fact ZQx corresponds to 
|(/)J)((/)J| — \(j)'^){<j>^\, and so both cat states wiU be re- 
tained by the triple quantum filter. The effect of revers- 
ing the network is then to convert this to the pseudo-pure 
state corresponding to 

|000)(000| - |100)(100| = 4 «) |00)(00|. (36) 

This is a pseudo-pure state of the last two spins, and in 
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general multiple quantum selection of cat states provides 
a convenient way of generating an n — 1 qubit pseudo- 
pure state in an n spin system. Furthermore, for some 
purposes states of the form given by equation B6| can be 



used as if they were n qubit pseudo-pure states |7! 



READOUT 

As described above there are two main methods for de- 
termining the final state of an NMR quantum computer: 
by examining the spectrum of the spins corresponding to 
the qubits of interest, and by examining the spectra of 
other neighbouring spins. These methods are simplest 
to describe when the quantum computer ends its com- 
putation with all the answer qubits in the eigenstates |0) 
and |1), rather than in superposition states or entangled 
states, as in this case a small number of measurements 
will provide all the information required |11| . A more 
thorough approach is to completely characterise the fi- 
nal state of the spin system by so called quantum state 
tomography [|2|; while the results can be interesting in 
small spin systems the effort required to perform tomog- 
raphy increases rapidly with the size of the spin system, 
and this approach is probably impractical for systems of 
more than three spins. 



Simple readout 

The simplest situation to consider is a one qubit 
NMR quantum computer which ends a calculation in the 
pseudo-pure state corresponding to |0) or |1). As dis- 
cussed above (equations ^ and |ll|) , these correspond to 
the NMR states Iz and — /^ respectively, and excitation 
with a 90° ly pulse will convert these to ±Ix- Thus the 
two states will give rise to absorption and emission lines 
in the NMR spectrum; this is hardly surprising as they 
correspond to excess population in the low energy and 
high energy spin states. It is, of course, necessary to 
obtain a reference signal against which the phase of the 
signal of interest can be determined, but this is easily 
achieved, either by using the NMR signal from a reference 
compound, or by acquiring a signal from the computer 
in a known state, |0) or |1). 

The situation is similar, but more complex, with larger 
spin systems. The NMR state corresponding to |00) is 
not just Iz + Sz, as might naively be expected; instead it 
is Iz + Sz + 2IzSz (see equation p^) . A general two qubit 
pseudo-pure eigenstate can be expressed similarly as 



\ab)(ab\ - i {{-iriz + i-lf Sz + {-ir®'2IzSz) 



spin, say /. The first approach is perhaps the most nat- 
ural approach in a homonuclear spin system, while the 
second method is more appropriate in a heteronuclear 
spin system. 

If both spins are excited, then the two population 
terms (Iz and Sz) are converted to single quantum co- 
herences, while the longitudinal two spin order is con- 
verted to unobservable double and zero quantum coher- 
ence. Thus the observable signal from a state of the form 
equation p^ is proportional to 



(-1)''/, + (-1)"^, 



(38) 



(37) 



This can be analysed in two ways: by exciting and ob- 
serving both spins, or by exciting and observing just one 



Clearly the desired information can be obtained from the 
phases (absorption or emission) of the NMR signals from 
the two spins. 

The situation is slightly more complicated if only one 
spin is observed: application of a 90° ly pulse to the state 
equation ^ gives 

((-1)"4 + (-1)"^, + {-ir®'2I,Sz) /2, (39) 

and the observable signal is proportional to 

i^ir {I, + i-lf2I,Sz) . (40) 

Thus only one of the two lines in the / spin doublet will 
be observed; which of the two lines this is depends on b, 
the state of spin 5", while the phase of the signal depends 
on a, the state of spin /, as before. 

Tomography 

Many NMR quantum computation experiments have 
used a readout scheme called quantum state tomography, 
and while this scheme is impractical for use with large 
spin systems it merits some explanation. The easiest ap- 
proach to readout is simply to determine the states of 
one or more critical qubits which contain the desired an- 
swer, while an alternative, far more thorough, approach 
is to characterise the complete density matrix describing 
the final state of the system [|l2[. This state tomogra- 
phy approach requires a large number of different mea- 
surements to fully characterise all the elements of den- 
sity matrix, and for large spin systems the complexity of 
this approach becomes prohibitive. For small systems, 
however, it provides detailed information not just on the 
result of the calculation, but also on any error terms. 

The density matrix describing a two spin system can 
itself be described using fifteen real numbers, correspond- 
ing to the amounts of the fifteen two-spin product opera- 
tors in the state (neglecting the identity matrix as usual). 
In a heteronuclear spin system it is possible to determine 
the values of four of these coefficients (the amounts of Ix, 
ly, 2Ij;Sz, and 2IySz) just by observing the I spin free in- 
duction decay, while four more can be determined by ob- 
serving spin S. The seven remaining coefficients can then 
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be determined in a minimum of two more experiments by 
exciting either I or S before observation. In general the 
spectrum of a single spin can provide at most 2" real 
numbers, while 4" ~ 1 numbers are required to charac- 
terise the spin system; thus a minimum of 2" separate 
experiments will be required. In practice the schemes 
actually used are substantially less efficient, greatly in- 
creasing the effort required for full tomography. For ex- 
ample, one tomographic analysis of a heteronuclear two 
qubit system involved nine separate experiments fTq. 



PRACTICALITIES 

Selective pulses 

Implementing these pulse sequences in a fully het- 
eronuclear spin system is straightforward, but in a 
homonuclear spin system complications arise from the 
need to perform selective excitation. The simplest ap- 
proach is the use of conventional selective pulses p5[ . 
These pulses can be simple Gaussian pulses incorporat- 
ing a phase ramp to allow off-resonance excitation, but it 
is probably better to use more subtle pulse shapes, such 
as members of the BURP family of pulses |Q . The soft 
pulses should excite all the lines in the target multiplet in 
an identical fashion, while leaving other lines completely 
untouched. In practice this is difficult to achieve in ^H 
systems, leading to the substantial errors clearly visible 
in many experiments. 

As pulse sequences implementing quantum logic gates 
can contain a large number of selective pulses separated 
by delays, it is necessary to address each spin in its own 
rotating frame. In homonuclear two-spin systems, how- 
ever, such as those used to implement two qubit NMR 
computers, it is possible to use a simpler approach. Sup- 
pose the centres of the two multiplets are separated by 
i/Hz; in this case the two frames will rotate with a rel- 
ative frequency i^. If the rotating frames were aligned 
at the beginning of the pulse sequence, they will come 
back into alignment at time intervals 1/v. As long as 
excitation and observation is performed stroboscopically 
it is possible to treat both nuclei as inhabiting the same 
rotating frame. Similarly, by choosing times such that 
the two rotating frames are 90 or 180° out of phase, it 
is possible to use variations on the simple "jump and re- 
turn" pulse sequence [[78| to perform selective excitation. 
This approach ||7^ can prove simpler than using selective 
pulses directly, but it cannot easily be used in systems 
with more than two spins of a given nuclear species. 



Composite Pulses 



Composite pulses [45 



play an important role in 



imental imperfections, such as pulse length errors and 
off-resonance effects, to be reduced. Such pulses could 
also prove useful in NMR quantum computers, acting to 
reduce systematic errors in quantum logic gates [ |81| . Un- 
fortunately most conventional composite pulse sequences 
are not appropriate for quantum computers as they only 
perform well for certain initial states, while pulse se- 
quences designed for quantum information processing 
must act as general rotors, that is they must perform 
well for any initial state. 

Composite pulses of this kind (sometimes called Class 
A composite pulses Q ) are rarely if ever needed for more 
conventional NMR experiments, and so have been rela- 
tively little studied. One important example is a compos- 
ite 90° pulse developed by Tycko |8^, |2|, which has re- 
cently been generalised to arbitrary rotation angles pu. 
These composite pulses give excellent compensation of 
off-resonance effects at small offset frequencies, such as 
those found for ^H nuclei, but are of no use for the much 
larger off-resonance frequencies typically found for ^^C. 

Fortunately when composite pulses are used for NMR 
quantum computation one great simplification can be 
made: it is only necessary that the pulse sequence per- 
form well over a small number of discrete frequency 
ranges, corresponding to the resonance frequencies of the 
nuclei used to implement qubits; it is not necessary to de- 
sign pulses which work well over a broad frequency range. 
In particular many NMR quantum computers use at most 
two spins of each nuclear species (see, for example, |75|), 
and it is convenient to place the RF frequency in the 
centre of the spectrum, so that the two spins have equal 
and opposite resonance offsets JT^l- Thus it is sufficient 
to tailor the composite pulse sequence to work well at 
these two frequencies, while the performance at all other 
frequencies can be completely ignored |83] . 



Abstract reference frames 

One technique which has proved extremely useful in 
the implementation of NMR quantum computers with 
more than two qubits is the use of abstract reference 
frames [ [72| . As it is necessary to address each spin in 
its own rotating frame of reference, it is possible to sim- 
ply rotate this frame to absorb the effects of z rotations, 
whether these arise from attempts to implement quan- 
tum gates (see equations ^ and ^), or the failure to 
fully refocus chemical shifts. 

For example, 905-^ rotations occur in the implemen- 
tations of many interesting gates. If need be these can 
be achieved either by periods of free precession, or by 
composite z-pulses. A simpler approach, however, is 
to achieve the same effect by rotating the RF reference 
frame, so that subsequent pulses are applied with appro- 
priate phase shifts. Thus, for example, the pulse sequence 



many NMR experiments, enabling the effects of exper- 90z 90x can be replaced by 90_j^ 90z: the phase of the RF 
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FIG. 10: A quantum circuit for the classical analysis of a bi- 
nary function /; this circuit is drawn for a quantum computer 
but is equivalent to that for a classical reversible computer 
(see figure m . 



pulse has been shifted, and the z-pulse has been delayed. 
Ideally it is possible to use this method to delay the z 
rotation to the very end of the pulse sequence, where it 
can be replaced by a rotation of the RF detection axis, 
or in many cases ignored all together. 



SIMPLE ALGORITHMS 

Now that we have seen all the elements necessary to 
implement quantum logic operations within NMR it is 
useful to see how they can be assembled to build small 
NMR quantum computers. Only two algorithms will be 
discussed in detail, both of which can be implemented us- 
ing two qubit computers, that is two spin systems. Brief 
reference will, however, be made to more complex sys- 
tems. 

Computers as small as these bear little immediate re- 
semblance to the computers in widespread use today: 
with only two qubits there is simply no available memory 
in which to store extraneous data or programs! Instead 
the program is built into the design of the NMR pulse se- 
quence used to implement the computation, and the two 
qubits are used to store the input data and the result of 
the computation, as well as forming the "CPU" of the 
system. 



Functions and phases 

Before discussing the algorithms themselves, it is useful 
to describe a trick widely used in quantum computation 
for converting the results of a function evaluation into a 
phase shift. This phase trick plays a central role in con- 
ventional implementations of many algorithms, but with 
NMR quantum computers it is often more appropriate 
to redesign the computer to implement the desired phase 
shifts directly. 

These simple demonstration algorithms are based on 
the analysis of one-bit binary functions, that is functions 
which take in one or more bits as input and return a 
single bit (that is, or 1) as output. These functions can 
be evaluated on reversible computers using /-controUed- 
NOT gates, as shown in figure o, with the result returned 
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TABLE I: The four possible binary functions mapping one bit 
to another; each function is conveniently labeled by the bit 
pattern in its truth table. 



as the value of an additional output bit, which begins 
the computation initialised to 0. An equivalent approach 
can be used with quantum computers, figure nG, but it 
is also possible to perform function evaluation with this 
"output" qubit set not to |0) but to the superposition 
(|0) - |l))/\/2. Since 
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(41) 



an /-controUed-NOT will perform the transformation 



k)(|o)-|i)) / 

V2 



(-1) 
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l^)(|o)-|i)) 
V2 



(42) 



and so the result of the function is returned as a phase. 
Note that the starting state of the ancilla qubit can be 
easily prepared from the state |1) by the application of a 
Hadamard gate (equation g) . 

This phase trick might seem pointless, indeed counter- 
productive, as it seems to return the result of the function 
as a global phase, and such global phases have no phys- 
ical meaning. As we shall see, however, the phase trick 
can be combined with quantum parallelism in a cunning 
and useful way. 



Deutsch's algorithm 

Deutsch's algorithm |2|, |^ g5| is concerned with the 
analysis of binary functions from one bit to one bit, that 
is functions which take in one bit as input and return 
another bit as output. Clearly there are four such func- 
tions, as shown in table |. These four functions can be 
divided into two groups: the two constant functions, for 
which f{x) is independent of x (/go and /n), and the two 
balanced functions, for which f{x) is zero for one value 
of X and one for the other (/oi and /lo). Equivalently, 
the functions can be classified according to the parity of 
the function, /(0)®/(l). 

Given some unknown function / (chosen from among 
these four functions), it is possible to determine which 
function it is by applying / to two inputs, and 1, using 
the circuit shown in figure Ru. This procedure also pro- 
vides enough information to determine the parity of /, 
and thus whether the function is constant or balanced. 
However knowing the parity of / corresponds to only one 
bit of information, and so it might be possible to answer 
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FIG. 11: A quantum circuit implementing Deutsch's algo- 
rithm to determine the parity of a binary function /. 



this question using only one evaluation of the function 
/. This cannot be achieved with a classical computer, 
but with a quantum computer the problem can be solved 
using Deutsch's algorithm. 

The basic idea behind Deutsch's algorithm is to com- 
bine the phase trick with quantum parallelism. Suppose 
that the /-controlled-NOT gate is applied with the input 
qubit in the state (|0) -I- \\))/y/2 and the ancilla qubit in 
the state (|0) — |l))/-\/2; then from equation ^ the result 
of the computation will be 



(_l)/(o)|0) + (_l)/(i)|i)W|0)-|l) 



V2 



V2 



= (-1) 



^./(O) / (|0) + (~1)/W^/W|1))(|0)-|1)) ^ 
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The ancilla qubit remains in (|0) — |l))/-\/2, while the "in- 
put" qubit now contains the state (|0) ± |l))/\/2, where 
the choice of plus or minus sign depends on /(O) © /(I). 
As relative phases in superpositions can be detected (for 
example, by applying a Hadamard gate as shown in equa- 
tions H and ^ , this allows the parity of / to be determined 
with only one function evaluation. 

A quantum circuit for Deutsch's algorithm is shown 
in figure O. The Hadamard gates act to interconvert 
eigenstates and superpositions, allowing both the phase 
trick and quantum parallelism to be implemented. Note 
that in this algorithm there is no input, and the result 
ends up in the first qubit, not the second qubit as occurs 
for traditional function evaluation. The second qubit is 
used simply as an ancilla to implement the phase trick. 
As discussed previously , it is not possible to program 
such a simple computer; instead the choice of function / 
is "hard wired" into the computer by the design of the 
/-controlled-NOT gate. 



NMR implementations 

Deutsch's algorithm was not only the first quantum 
algorithm to be described: it was also the first algo- 
rithm to be implemented on an NMR quantum computer, 
first using a homonuclear system (cytosine) and then a 
heteronuclear system (^'^C-labeled chloroform). These 
two implementations will be described below, while more 
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FIG. 12: Modified quantum circuit for the classical analysis of 
/(O) on an NMR quantum computer; Uf is a propagator cor- 
responding to the /-controlled-NOT gate in the conventional 
circuit (figure hol). Function evaluation is followed by 90y 
pulses to excite the NMR spectrum. Clearly /(I) can be ob- 
tained in a very similar fashion. 
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FIG. 13: Modified quantum circuit for the implementation 
of the Deutsch algorithm on an NMR quantum computer. 
Hadamard gates have been replaced by pseudo-Hadamard and 
inverse pseudo-Hadamard gates, that is 905;^ pulses. The final 
90y excitation pulses cancel out the 901^ pulses, and thus all 
four pulses can be omitted. 



modern implementations, including some extensions and 
simplifications, are described in the next section. 

The cytosine system ]1^ used the two ^H nuclei re- 
maining on a cytosine molecule when dissolved in D2O 
(figure 0). The two ^H multiplets are separated by 763 Hz 
(at a ^H frequency of 500 MHz), with a J-coupling of 
7.2 Hz. Selective excitation was achieved using Gaussian 
shaped soft pulses incorporating a phase ramp, and the 
lengths of the soft pulses were chosen as multiples of the 
inverse of the frequency separation of the two resonances, 
so that the unexcited spin experienced no net rotation 
during a selective pulse. 

Both classical function analysis and Deutsch's algo- 
rithm were implemented, using the modified quantum 
circuits shown in figures ^ and |l^. In these figures the 
/-controlled-NOT gates have been written as general two 
qubit propagators, C//, Hadamard gates have been re- 
placed by 90j_y pulses, and final 90° pulses have been 
added at the end of each sequence to convert eigenstates 
to observable magnetisation; for simplicity pulses which 
simply act to cancel one another are omitted. The ini- 
tial pseudo-pure states were prepared using field gradient 
techniques (spatial averaging), and the final results were 
determined by examining the phase (ix) of the NMR 
signals from the two spins. 

The heart of these computations is contained in the 
implementation of the propagators Uf . Each propagator 
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corresponds to flipping the state of the second spin as 
follows: C/oOj never flip the second spin; Uqi, flip the 
second spin when the first spin is in state |1); Uio, flip 
the second spin when the first spin is in state |0); Un, 
always flip the second spin. The first and last cases are 
particularly simple, as t/oo corresponds to doing nothing, 
while Uii is just a selective 180° pulse (a not gate) on 
the second spin. The second propagator is a controUed- 
NOT gate, while the third case is a reverse controlled-NOT 
gate, so that a not gate is applied to the target when the 
control spin is in state |0). These last two gates can be 
implemented as described above. The sequences actually 
used were 



1 
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90Sy - 77 - 180, - 77 - 180^90Iy90I,90.y90S±^ 
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AJ 
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where pulses not marked as either I or S were applied 
to both nuclei. The phase of the final pulse distinguishes 
Uqi (for which the final pulse was S^x) from Uio (for 
which it was S-x)- In retrospect it is clear that these se- 
quences are unnecessarily complicated; substantial sim- 
plifications could have been achieved by combing pulses 
and by absorbing phase shifts into abstract reference 
frames. 

The second implementation of Deutsch's algorithm 
used the heteronuclear two spin system provided by ^"Re- 
labeled chloroform in solution in deuterated acetone. The 
initial pseudo-pure state was prepared by temporal av- 
eraging, although results were also shown for computa- 
tions beginning in the thermal equilibrium state, while 
the results of the computation were determined both by 
direct observation of the ^H spectrum, and by full quan- 
tum state tomography, which allows errors to be studied. 
Only the Deutsch algorithm was demonstrated, with no 
results shown for classical computations. 

The pulse sequences used to implement the Uqi and 
Uio were similar to those used for cytosine; slight dif- 
ferences can be traced to the heteronuclear nature of the 
spin system and the consequent ability to place both spins 
exactly on resonance. The sequences used to implement 
the two balanced functions, f/oo and Uu were more com- 
plicated than for cytosine, as periods of free precession 
were included so that each of the four propagators was 
applied over approximately the same time period (about 
1/2 J). This means that the effects of relaxation (domi- 
nated by the relatively short T2 of the ^^C nucleus) were 
similar in all cases. 



Extensions and simplifications 

In addition to these two early examples several more 
implementations of Deutsch's algorithm, and its more 
general cousin the Deutsch-Jozsa algorithm [E3, have 



been published. From among these a few particularly 
interesting examples will be described in more detail. 

The Deutsch- Jozsa algorithm is a generalisation of 
Deutsch's algorithm which considers binary functions 
with any number of input bits. Clearly such functions 
need not be either constant (that is, give the same out- 
put for all input values) or balanced (output for half 
the possible inputs and 1 for the remainder); for exam- 
ple a binary function with two input bits could return 
for one of the four input values and 1 for the other 
three. Suppose, however, that it is guaranteed that some 
(otherwise unknown) function / is either constant or bal- 
anced (such theoretically convenient if apparently arbi- 
trary guarantees are usually referred to as promises), and 
it is necessary to determine to which of these two cate- 
gories the function belongs. To solve this problem by 
classical means would in the worst case require the evalu- 
ation of / for just over half its inputs (if the function uses 
n input bits it may be necessary to evaluate the function 
over 2"~^ -I- 1 inputs), while even in the best case it would 
be necessary to evaluate the function at least twice. By 
contrast the Deutsch-Jozsa algorithm can always distin- 
guish between constant and balanced functions with a 
single function evaluation. 

The first experimental implementation of the Deutsch- 
Jozsa algorithm |Q, which used functions with two in- 
put bits, is also notable as the first example of the use of 
transition selective pulses. The three qubit system cho- 
sen was the homonuclear spin system made up by the 
three ^H nuclei in 2,3-dibromopropanoic acid. For sim- 
plicity pseudo-pure states were not prepared; instead the 
algorithm was simply applied to the spin system in its 
thermal state. Function evaluation for the six possible 
balanced functions was accomplished by the simultane- 
ous application of single transition selective 180° pulses 
to two of the four components of the low field multiplet 
(which exhibits the largest separations between the four 
components), and the result of the computation was de- 
termined by observation of NMR signal intensities in the 
two high field multiplets. This simple approach worked 
remarkably well; the deviations from ideal behaviour seen 
in experimental spectra were largely ascribed to the ef- 
fects of strong coupling. 

The Deutsch-Jozsa algorithm has also been applied to 
larger spin systems, most notably a largely heteronuclear 
5 spin system (containing single ^H, ^^N and ^^F nuclei, 
together with two ^^C nuclei) derived from glycine [ [75[ ; 
in this case multiple pulse techniques were used to im- 
plement controlled gates. This system permits functions 
with 4 input bits to be studied, but only one constant 
and one (particularly simple) balanced function were ac- 
tually implemented, out of a possible total of two con- 
stant and 12870 balanced functions. Similarly while an 
approach related to temporal averaging was used for ini- 
tialisation, it was used to construct an initial state which 
gave the same signal as a pseudo-pure state, rather than 
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the pseudo-pure state itself. Thus this can only be con- 
sidered as a partial implementation. 

Another important variation on the Deutsch-Jozsa al- 
gorithm is to remove the ancilla qubit which is normally 
used to implement the phase trick, converting function 
values into phase shifts. This indirect approach is not 
necessary for NMR implementations, as it is possible to 
implement controlled phase-shift gates directly. Indeed it 
is usually simpler to do this; in particular the propagator 
for both constant functions is reduced to "do nothing" . 
This approach, sometimes called the refined Deutsch- 
Jozsa algorithm [|^, has the advantage that one fewer 
qubit is required to implement a given algorithm. It does 
however have one minor disadvantage for demonstration 
algorithms, as such systems are intrinsically quantum 
mechanical and cannot be used to implement classical 
function analysis. 

This simplified approach has been used with a three 
spin system (the three ^'^C nuclei in labeled alanine) 
to implement the Deutsch-Jozsa algorithm for functions 
with three bit inputs |pq|. In this case there are 70 bal- 
anced functions, from which ten representative functions 
were chosen. At the other extreme this technique can 
also be used to implement the refined Deutsch algorithm 
using a single qubit! In this case the algorithm is so 
simple as to be almost trivial: the Uf propagator is a 
180° pulse for the two balanced functions, while for the 
two constant functions it is as usual "do nothing" . Since 
the z-pulse can be absorbed into the RF reference frame, 
the pulse sequence can in principle be reduced to a 90° 
pulse followed by observation along ±x. Such a simple 
experiment seems hardly worth performing, but for com- 
pleteness it has been demonstrated as one member of a 
set of experiments using the isolated one spin and two 
spin ^H systems is 5-nitro-2-furaldehyde jSq]. 



Grover's quantum search 

Grover's algorithm |8^, Q is designed to speed up 
searches comparable to searching for a needle in a 
haystack. More mathematically it concerns the analysis 
of binary functions which map a large number of bits to a 
single output bit, where the task is to determine an input 
for which the value of the function is 1. If the function 
has many inputs for which its value is 1 (that is, if the 
haystack contains a substantial number of needles), one 
of these inputs can be easily located by trial and error, 
but if there is only one suitable input among a large num- 
ber of unsuitable inputs (one needle in a large haystack), 
locating this single input is obviously a difhcult process. 

Suppose that the function / has inputs described by 
n bits, so that there are N — 2" possible inputs, and 
that f — 1 for only one of these inputs. The only general 
way to locate this input is to evaluate / over some trial 
inputs, and look for a value / = 1 (a satisfying input). 



A lucky guess would permit this input to be located in 
one try, but this can hardly be relied on; on average a 
random search (the best classical algorithm) would re- 
quire about N/2 trial evaluations, or A^ — 1 evaluations 
in the worst case. The situation is similar if there are 
k inputs which satisfy the function; in this case about 
N/k trials will be required. By contrast, Grover's algo- 
rithm permits a satisfying value to be located with only 
0{^J N/k) evaluations. This increase in computational 
efficiency is less impressive than the exponential increase 
seen for Shor's quantum factoring algorithm H, but is 
still quite substantial; furthermore the algorithm is quite 
general, with a range of potential applications. 

Early NMR implementations [^^ g9| ^ concentrated 
on the case n = 2, so that there are 4 possible inputs, with 
only a single satisfying input. For this case the operation 
of the algorithm is fairly simple to explain, and more 
complicated cases can be understood by analogy. The 
algorithm involves two steps: evaluation of the function 
over all possible inputs, followed by a selection process to 
pick out the desired result. As an example I will assume 
that /(Ol) = 1 while f{x) = for all other inputs x. 

The algorithm begins with one quantum register (that 
is a group of qubits) in a uniform superposition of the 
four possible inputs, so that its state is 



(|00) + |01) + |10) + |ll))/2. 



(45) 



An attempt to read out the value of this register will re- 
turn one of the four possible inputs at random. A prop- 
agator implementing the function / is then applied, so 
that / is evaluated over all 4 inputs; the propagator is 
set (either directly or indirectly by means of the phase 
trick) to return the value of / as a phase shift, that is 



Uf 



\x)^{-iy^'^->\x), 
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so that at the end of the calculation the quantum register 
is in the state 



(|00)-|01) + |10) + |ll))/2. 



(47) 



The desired satisfying input has now in some sense been 
identified, as it bears the unique mark of a negative 
phase; this is not, however, of any immediate use, as 
an attempt to analyse this state will still return one of 
the four contributing inputs at random. It is, therefore, 
necessary to apply some propagator which converts the 
phase difference into an amplitude difference. 

This process might seem simple, but it is in fact quite 
tricky, as any such propagator must correspond to a log- 
ically reversible unitary operation. There is, however, a 
solution: inversion around the average. This slightly pe- 
culiar operation takes in a superposition and reflects the 
amplitude of each component around the average ampli- 
tude of all the components. In the example, the individ- 
ual amplitudes are ±i, and the average amplitude of the 
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FIG. 14: A quantum circuit for the implementation of 
Grover's quantum search algorithm on a two qubit computer. 
Boxes marked H are Hadamard gates. The first two qubit 
gate Uf^^ corresponds to evaluation of the function fab, re- 
placing an eigenstate \ij) by —\ij) ii i — a and j = 6, while 
Ufgg simply replaces |00) by — 100). 



four components is j; reflecting ^ around j gives 0, while 
reflecting — ^ gives 1. Thus this operation acts to concen- 
trate all the amplitude on one member of the superposi- 
tion, giving a final state of just |01). Surprisingly this ap- 
parently complex operation can be performed using only 
Hadamard gates and a controlled phase-shift gate which 
negates |00) while leaving all other states alone. A quan- 
tum network for implementing this algorithm is shown 
in figure R% this network assumes that the function is 
evaluated using propagators which apply the necessary 
phase shifts directly; this can of course be achieved using 
an ancilla qubit if desired, but direct application of phase 
shifts has been the universal practice in NMR implemen- 
tations of this algorithm. 



final pulse sequences for Uab were 



^ ~ 90_j,905±,90/±, -^- 90_y90_, (48) 

where, as before, pulses not marked as either I oi S were 
applied to both nuclei. The choice of ± signs on the 
second and third pulses determined which of the three 
functions / is implemented. Results were shown only for 
Uii, in which case both signs are positive, but experi- 
ments were performed for all four functions. 

The second NMR implementation of Grover's quantum 
search [g9[ was based on the homo nuclear ^H spin sys- 
tem in cytosine. In this case spatial averaging was used 
to prepare the initial pseudo-pure state, and the result 
was analysed by excitation and detection of the two "'^H 
signals; this provides a particularly simple and immedi- 
ate readout scheme. Quantum logic gates were imple- 
mented using multiple pulse sequences, with soft pulses 
used to perform selective excitation and spin echoes used 
to refocus chemical shifts; sequential pairs of pulses were 
combined within the individual propagators f//^^, but no 
attempt was made at global simplification of the pulse 
sequence. For these reasons the homonuclear implemen- 
tation produced relatively poor results, although the cos- 
metic appearance of the spectra was greatly improved by 
the application of a magnetic field crush gradient between 
the end of the quantum circuit and the application of a 
final 90° pulse prior to detection. 



Extensions 



NMR implementations 

The first implementation of Grover's search |7^ was 
performed using a heteronuclear NMR quantum com- 
puter based on chloroform. This used temporal averag- 
ing to prepare the initial pseudo-pure state, and quantum 
state tomography to characterise the final result. Quan- 
tum logic gates were implemented using multiple pulse 
sequences as discussed below. 

As this implementation used a heteronuclear spin sys- 
tem, both nuclei were placed on resonance in their respec- 
tive rotating frames; thus it was not necessary to refocus 
chemical shifts, and periods of free precession correspond 
to evolution under the spin-spin coupling. The four de- 
sired controlled phase shift gates were achieved (up to 
an irrelevant global phase) by combining this with ±z 
rotations on the two spins; these were explicitly imple- 
mented using composite z-pulses. The Hadamard gates 
were implemented using the two pulse sequence 90_yl80:E 
described in section . Finally the elements of the pulse se- 
quence with the exception of the initial pair of Hadamard 
gates were assembled together to give a single propagator, 
Uab, and sequential pairs of pulses were combined where 
ever possible to give simpler pulse sequences. Thus the 



Grover's quantum search algorithm is not, of course, 
limited to two qubit implementations, but can be used 
to search over a space described by any number of input 
qubits, n, in which case there are N = 2^ possible in- 
puts. For n > 2 the behaviour of the algorithm is similar 
to but more complex than that described above; a sin- 
gle application of Uab (that is, function evaluation and 
inversion around the average) acts to concentrate the in- 
tensity of the superposition on the satisfying state, but 
does not simply produce this state. Instead it is nec- 
essary to apply these two operations repeatedly, driving 
the register towards the desired state. The intensity of 
the desired state oscillates with a frequency inversely pro- 
portional to VN , and so after 0('\fN) applications of Uab 
the intensity will be largely on the satisfying state; mea- 
surement of the register will then return this state with 
high probability. Further application of Uab will then 
drive the register away from the desired state, and so it 
is important to choose the correct number of repetitions. 
This oscillatory behaviour was in fact demonstrated for 
the two qubit case in the first NMR implementation ^j\ , 
where Uab was applied between zero and seven times. 
More recently Grover's algorithm has been implemented 
with three qubits, searching over eight possible inputs. 
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with up to 28 repetitions of the propagator [£1| ; this im- 
plementation used the ^H-^^F-^'^C spin system in ^"^C 
labeled CHFBra. 

Another variant on Grover's algorithm occurs when 
there is more than one satisfying input; in this case the 
number of such inputs is usually called k. The algorithm 
is almost identical to the simple case when A; = 1, except 
that it is only necessary of use 0{y/N/k) repetitions to 
drive the quantum register into an equally weighted su- 
perposition of the k satisfying inputs. A measurement on 
this register will then cause the superposition to collapse 
into one of its constituent values, and so one of the k sat- 
isfying inputs can be returned at random. Clearly this 
requires either that the value of k be known beforehand, 
or that it be determined; fortunately k can be readily 
estimated using an extension of Grover's search called 
quantum counting [|9[ || || ||. 

The basic idea behind quantum counting is that in ad- 
dition to driving the quantum register towards the sat- 
isfying values application of the Grover propagators also 
results in a phase shift which depends on the value of 
k/N. In a conventional implementation of Grover's al- 
gorithm this phase shift is a global phase shift, and so 
cannot be detected. Quantum counting, however, uses 
an approach similar to Deutsch's algorithm to convert 
this phase shift into a relative phase shift which can be 
measured. This algorithm requires one additional qubit, 
and so when it was implemented on the cytosine system 
p9[ it was performed using one bit functions, for which 
either zero, one, or two inputs satisfy the function. 



QUANTUM PHENOMENA 

In addition to implementing quantum computations, 
NMR techniques have also been used to implement other 
more general aspects of quantum information processing, 
including demonstrations of some quantum phenomena. 
A few of the more important examples are described be- 
low. 



Entangled states 

Entangled states [Tg] are states of quantum mechani- 
cal systems which cannot adequately be characterised by 
describing the states of their component subsystems; in- 
stead the properties of such states are properties of the 
system as a whole. An important example is provided by 
the four Bell states: 



|00)±|11) 
V2 



^' 



|io)±|oi) 

V2 ■ 



(49) 



quantum theory, and also form the essential information 
processing resource for quantum communication tech- 
niques. For this reason there is significant interest in 
techniques for generating entangled systems, especially 
multiple particle entangled systems such as cat states. 

The quantum network for generating cat states is well 
known, and the three qubit version is shown in figure 
^. This was used early on to generate Bell states and 
three qubit cat states (GHZ states) Q, and more re- 
cently has been used to prepare seven qubit cat states 
Q. There are, however, several reasons for question- 
ing the true significance of these results, and they have 
certainly not generated as much excitement as similar re- 
sults with smaller numbers of qubits in other quantum 
technologies 
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The first reason for skepticism is the close relationship 
between cat states and multiple quantum coherences. As 
described above (Section ) nQx multiple quantum coher- 
ence corresponds to a mixture of \4>n){4'n\ ^'^^ \4>n){4>n\'i 
thus the generation of high order multiple quantum co- 
herence is nearly equivalent to the production of cat 
states. Seen in this light, seven quantum coherence is 
not particularly impressive: solid state NMR techniques 
have been used to generate coherence orders above one 
hundred ^. 

A second difficulty with NMR cat states arises from 
the use of pseudo-pure states. As discussed in Section 
, the fact that NMR density matrices are always highly 
mixed, with nearly equal populations in all spin states, 
appears to mean that they cannot, strictly speaking, ex- 
hibit entanglement. Thus NMR cat states might be more 
properly described as pseudo- entangled states. 

Finally even if truly entangled NMR states were to be 
produced there are are serious limitations on the use of 
NMR to investigate quantum mechanics. The ensemble 
nature of NMR measurements complicates the investi- 
gation of deviations from classical behaviour, while the 
short distances over which NMR entanglement can be 
produced (normally confined to molecular dimensions) 
compares unfavourably with the distance achievable with 
entangled photons (hundreds of metres). 



Quantum teleportation 



(The cat states introduced in Section are an obvious gen- 
eralisation of <f) to systems of more than two qubits.) 
Such states play a central role in experimental tests of 



Quantum teleportation |^, ^ |100| is a particu- 
larly intriguing example of quantum communication; in 
essence it involves the transfer of an unknown quantum 
state from one quantum particle to another, without any 
attempt to characterise the state. The technique relies on 
the peculiar, apparently non-local, correlations inherent 
in entangled states, such as Bell states. Note that quan- 
tum teleportation does not permit the direct transfer of a 
quantum system into empty space: a suitable target par- 
ticle (one half of a Bell state) must be provided at the 
destination to receive the quantum information. Thus it 
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is the state of the particle which is teleported, and not 
the particle itself. 

Quantum teleportation has been implemented on a 
three spin NMR quantum computer |101] using two "'^^C 
nuclei and a single ^H nucleus in ^'^C labeled trichloro- 
ethene. The process can be summarised by 

IV'a) (lObOe) + \lblc)) — ^ (lOaOfc) + |lal&» \A) (50) 

where 1-0) indicates an arbitrary quantum state, and the 
subscripts a, b, c simply label the three qubits. Although 
this implementation raises some interesting issues, argu- 
ments similar to those used above for entangled states can 
be applied to NMR teleportation, and once again optical 
implementations |100[ are rather more convincing. 



Error correction 

Any computing technology is ultimately based on some 
physical device, and such devices are inevitably error 
prone. There are two main methods by which the ef- 
fects of these random errors can be reduced: stabilization 
techniques which act to cancel out the effects of small er- 
rors, and error correction techniques which detect, char- 
acterise and finally fix the results of larger errors. 

Stabilization against small errors is an inherent feature 
of digital information processing. In any digital system 
information is stored as ones and zeroes, which are ul- 
timately represented as two states of a physical system, 
such as high and low voltages. Small fluctuations away 
from the two ideal voltages (noise) do not matter, as long 
as it is always clear whether the voltage is high or low. In 
some cases this passive insensitivity to noise can be fur- 
ther enhanced by active stabilization, which acts to con- 
tinuously drive the system towards the nearer of the two 
ideal states. Because of this intrinsic stabilization digital 
information processing devices are effectively invulnera- 
ble to the effects of noise, as long as the noise signals 
remain below some critical threshold. If the noise rises 
above this threshold, however, stabilization is no longer 
effective, and it is necessary to resort to error correction. 

Error correction techniques are even more important 
for quantum information processing, as simple stabiliza- 
tion techniques are ruled out. Unlike bits qubits are not 
confined to two states, but can also exist in superposi- 
tions of these states, and any stabilization scheme which 
drives a qubit back towards the two eigenstates will of 
course destroy these vital superpositions. For some time 
it was believed that the nature of superposition states 
would also render error correction schemes impractical, 
but happily this is not in fact the case. 

Classical error correction schemes [15| are most sim- 
ply described in terms of the transmission of information 
along a noisy channel, which has the effect of flipping 
bits (that is, changing them from to 1 and vice versa) 



at random, with an error probability e. Quantum error 
correction schemes are similar except that as well as cor- 
recting qubit flip errors (that is errors which take a qubit 
from |0) to |1)) it is also necessary to correct phase errors 
(which interconvert, for example, |0) -I- |1) and |0) — |1)). 
For simplicity, however, I will only describe how classi- 
cal schemes work. All such schemes (including quantum 
schemes [102, 103, 104, 105|) use multiple copies of each 
bit (or qubit); the redundancy provided by these ancillas 
allows errors to be detected and corrected. 

The simplest classical scheme is triplet coding, in which 
a bit is encoded using three repetitions (so that is en- 
coded as 000, while 1 is encoded as 111) and decoded by 
taking a majority vote (so that 000, 001, 010 and 100 
all decode as 0, while 111, 110, 101 and Oil all decode 
as 1). This scheme is robust against random errors in 
any one of the three bits. Of course if there are errors in 
two of the bits then the message will still be corrupted; 
however the chance of two errors occurring is e^(l — e), 
and as long as e is small this possibility can be neglected. 
If the level of errors is too high then triplet coding is no 
longer effective, but more robust schemes (involving even 
greater redundancy) can be used. 

Triplet coding and other error correction schemes 
might seem very different from stabilization schemes, but 
in fact the basic ideas are quite similar as shown in figure 
|l5| . In effect the code divides the eight possible settings 
of a three bit system into two subspaces, just as voltages 
can be divided into high and low. Note that if a simpler 
doublet coding scheme (in which is represented as 00 
and 1 as 11) is used it is still possible to detect single bit 
errors, but not to correct them, as the two states which 
can occur after a bit flip error (01 and 10) lie equally close 
to the two ideal states. In communication (as opposed 
to data storage) schemes, however, it may be sufficient 
to detect errors, as the erroneous bits can then be sent 
again. 

Some simple quantum error detection and correction 
protocols have been implemented on NMR quantum com- 



puters 1 106, 107 1 . Full quantum error correction requires 
at least five qubits to encode a single state, but simpler 
schemes exist which use only three qubits; these sim- 
plified schemes can only correct phase errors or bit-flip 
errors, but not both. Phase errors occur as a result of 
spin-spin relaxation, while bit-flip errors correspond to 
spin-lattice relaxation, and so in many NMR systems 
phase errors will dominate. Furthermore, phase errors 
only occur in quantum computers, as they have no clas- 
sical analogue. For these reasons early studies on NMR 
error correction have concentrated on three qubit phase- 
correcting codes. 

It should be noted that these NMR experiments are 
demonstrations of the principle of error correction, rather 
than practical implementations of error correcting codes. 
In order to effectively suppress errors in a quantum com- 
putation it is necessary to apply the error correction pro- 
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FIG. 15: The relationship between the triplet error correc- 
tion code and active stabilization. The eight possible states 
of three bits can be placed on the eight corners of a cube, 
where the sides of the cube connect states which differ by a 
single bit flip and the two ideal states (000 and 111) lie at 
opposite corners. Decoding a triplet state by majority vote is 
equivalent to moving the state to the nearer of the two ideal 
corners. 



tocol repeatedly; this in turn requires that the ancilla 
qubits be maintained in the correct state. This is most 
simply achieved by initiaUsing them to |0) before each 
correction round. Unfortunately the NMR techniques 
described in Section allow qubits to be initialised only 
once, at the start of the calculation; they cannot be re- 
peatedly reinitialised. This appears to rule out current 
NMR implementations as practical technologies for quan- 
tum computation |B8|. 



NMR AND ENTANGLEMENT 

Finally I will return to the question, briefly discussed 
in Section , of whether NMR quantum computers are in 
fact true quantum computers at all. Much of the op- 
position to NMR as a quantum computing technology 
stems from the formidable difficulties [Q in scaling up 
the current small systems to computers with a reasonable 
number of qubits. A particularly common observation is 
that the use of pseudo-pure states is exponentially inef- 
ficient IQ ^ , in that the amount of pseudo-pure state 
which can be extracted from an NMR system at ther- 
mal equilibrium falls off exponentially with the number 
of spins in the system. Of course this problem can be 
overcome by using an exponentially large sample, but this 
approach would remove any increase in computational ef- 
ficiency supposedly arising from quantum mechanical ef- 
fects: there are a wide range of classical techniques (such 
as DNA computing ]10^) which allow exponential gains 



in computing power to be obtained from exponentially 
large samples. 

This problem is not in principle unique to NMR; it will 
occur in any potential quantum computing technology 
which works in the high temperature limit [Q . It can in 
principle be overcome by working at sufhcicntly low tem- 
peratures, or by using some other initialisation technique 
to produce a non-Boltzmann population distribution, al- 
though the technical problems involved are substantial 
p8[ . However NMR is the only technology among those 
currently under investigation which falls into this cate- 
gory. 

In addition to the obvious technological issues raised by 
this exponential efficiency there are also some more fun- 
damental concerns. It has long been suspected that the 
non-classical power of quantum computation is closely 
linked to the existence of entangled states during quan- 
tum computations [109|. Although this belief has never 
actually been proved, and some recent theoretical re- 
sults have suggested that it may not be entirely correct 
1 110 1, it is clear that some important algorithms such as 
Shor's quantum factoring algorithm do require the gen- 
eration of entangled states |111|. It can be shown that 
the pseudo-entangled states generated in NMR quantum 
computations do not actually fulfill the mathematical re- 



quirements for true entanglement [112|, casting doubt on 
their ability to exhibit true quantum phenomena. As we 
shall see, however, this concern may not be entirely well 
founded. 



Quantifying entanglement 

Although the concept of entanglement is relatively easy 
to explain, actually quantifying the amount of entangle- 
ment in any given system is a surprisingly difficult task. 
For a system of two qubits in a pure state the problem 
is relatively straightforward, and the four Bell states (0^ 
and tp^, see Section ) form a basis set describing the pos- 
sible maximally entangled states. With larger systems 
the problem is much more difficult, as different defini- 
tions of entanglement lead to quite different conclusions. 

Difficulties can also arise when considering mixed 
states, such as those observed in NMR experiments. 
These difficulties occur because there is no unique way 
to break down a given mixed state into a mixture of pure 
states. To see this consider the maximally mixed state, 
which for a two qubit system has the form 
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This is perhaps most easily described as an equally pop- 
ulated mixture of the four eigenstates, but this is by no 
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means the only possible description: an equally weighted 
mixture of the four Bell states will have exactly the same 
form! It might be claimed that this alternative decom- 
position is unnatural, but there are no real grounds for 
such a statement, as the choice of basis set is entirely 
arbitrary, and this approach is sometimes referred to as 
the preferred ensemble fallacy. 

Similar difficulties arise when considering the pseudo- 
entangled states formed from pseudo-pure states, such 
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(mix ed states of this kind were first considered by Werner 
1 113 1, and thus are often referred to as Werner states). It 
is tempting to argue that this mixed state is a mixture 
of the maximally mixed state together with a fraction e 
of an entangled state, but as discussed above there is no 
particular reason to choose this description. The amount 
of apparent entanglement in the state could be increased 
by decomposing the maximally mixed state as a mix- 
ture of entangled states, but it is also possible to choose 
decompositions which reduce the apparent contribution 
from entanglement. 

While it is not possible to say how much entanglement 
is in any particular mixed state it is possible to determine 
the minimum contribution from entangled states which 
must be present in the state (in a pure entangled state 
this minimum fraction is, of course, one). If this mini- 
mum quantity is greater than zero then it is reasonable 
to say that the mixed state does contain some entangle- 
ment; if, however, the minimum amount is zero then it is 
possible to describe the mixed state using only product 
states (and thus without invoking entanglement), and the 
state is said to be separable. In particular it can be shown 
1 114 1 that pseudo-entangled states are in fact separable if 
e < 1/3. An explicit separable decomposition of equation 
B2 with e = 1/9 is given in appendix . 

With larger systems the problem is more complicated, 
but two important results are known |112]. Given any 
state I'i/') of an n-qubit system, a mixture made by mixing 
a fraction e of this state into the maximally mixed state, 
1/2", can always be shown to be separable for sufficiently 
small values of e, such that 



e < 



1 + 22"-i 



(53) 



It can also be shown that non-separable states do exist 
for sufficiently large values of e, such that 



e > 
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By comparison the values of e obtainable with NMR 
quantum computers working within the high tempera- 
ture limit dl, |6[ 13 are given by 



n 
2" 



(55) 



which lies between the two bounds given in equations |53| 
and ^. Using realistic parameters it can be calculated 
that the states used in NMR quantum computations are 
always separable for systems with less than about 13 
qubits and may (or may not) become entangled beyond 
this point. Since the systems used so far have involved no 
more that seven qubits, all NMR quantum computations 
to date have involved purely separable states. 



NMR and quantum mechanics 

The observation that NMR quantum computers have 
so far only used separable states has led some authors to 
suggest that they are not true quantum computers at all! 
When assessing claims of this kind it must be remem- 
bered that quantum is being used here in its technical 
sense of provably non- classical. Consider, for example, a 
set of NMR quantum computers which are identical ex- 
cept for having different values of e, with some of them 
lying above the entanglement limit discussed previously, 
and the remainder lying below this limit. It seems very 
strange to claim that two groups of computers are funda- 
mentally different in character, with the first group being 
quantum mechanical and the second group classical, but 
it is more reasonable to suggest that only the computers 
in the first group are capable of exhibiting convincingly 
non-classical behaviour. 

Even this claim, however, may be too strong. It seems 
highly unlikely that it is the mere presence of entangle- 
ment which leads to non-classical efhciency; rather it is 
the ability to interconvert a wide range of entangled and 
non-entangled states. Thus in order to claim that NMR 
quantum computing experiments are classical it is not 
sufficient to show that they involve only classical states; 
instead it is necessary to show that the processes which 
connect these states can themselves be described classi- 



cally. To date attempts to achieve this have failed [115|, 
and it is not clear that such a model is possible. 

In an unrelated approach, some authors have at- 
tempted to draw a distinction between the density matrix 
(which is a description of the state of an NMR system) 
and the state itself, although it is tricky to draw this 
distinction in an entirely convincing fashion. It is true, 
however, that the density matrix approach is only an ap- 
proximate description of an NMR system, and that any 
conclusions based on this approximation are to some ex- 
tent open to suspicion. 
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CONCLUSIONS 



Acknowledgements 



When assessing NMR quantum computation it is im- 
portant to take a balanced view, avoiding both excessive 
excitement at the apparently impressive results achieved 
so far and undue despair at the limitations that have 
been identified. NMR quantum computation has been 
the subject of a great deal of skeptical scrutiny, probably 
more than any other approach. In part this is a result of 
the great success of NMR as a technique for quantum in- 
formation processing; furthermore, the highly developed 
nature of NMR experiments, in comparison with many 
other putative quantum technologies, means that the lim- 
its of the technique are well known and understood. 

On the positive side, NMR is far ahead of any compet- 
ing technology in the implementation of quantum com- 
putations and other forms of quantum information pro- 
cessing. Although some basic elements have been im- 
plemented using othe r tec hnologies, such as the ion trap 
controlled-NOT gate 116 , NMR remains the only tech- 
nology capable of implementing a complete quantum al- 
gorithm. Progress from two qubit devices [|[ |[ |l^, |ll|, |l^ 
to systems with seven qubits Q has been extremely 
rapid, and there is every reason to believe that more 
progress will soon be made. 

Against this it must be pointed out that most re- 
searchers believe that the current designs for NMR quan- 
tum computers cannot be extended very much farther: 
while there is some disagreement as to which technical 
difficulty will actually stop further progress it is widely 
agreed that it will be difficult to progress beyond 10- 
20 qubits. While current demonstration systems are un- 
doubtedly interesting they are far too small to be used 
to tackle problems beyond the range of current classi- 
cal computers. Similarly, although there are important 
applications of quantum info rmation processing, such as 
quantum cryptography | 117 ], which require only devices 
with small numbers of qubits, all such applications lie in 
the field of quantum communication where NMR meth- 
ods appear completely unsuitable. 



The discovery [112| that current NMR implementa- 
tions of quantum computation do not seem to involve 
entanglement might appear a serious blow, but its impli- 
cations should not be overstated. This result does not 
mean that NMR quantum computers are not true quan- 
tum computers, although it does appear to mean that 
they cannot be used to achieve non-classical eflficiencies. 
However it has long been known Efl, M that the ex- 
ponential inefficiency of pseudo-pure state preparation 
means that current implementations are unlikely to ex- 
hibit true efficiency gains. It seems likely that in the 
next few years it will be possible to use para-hydrogen 
techniques |^, |3|, |^ to build two qubit NMR quantum 
computers above the entanglement threshold, but it will 
be tricky to extend this approach to larger systems. 
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AN EXPLICITLY SEPARABLE 

DECOMPOSITION OF A PSEUDO-ENTANGLED 

STATE 

While it can be shown that a two qubit pseudo- 



entangled state with e < 1/3 is in fact separable |114|, the 
argument used to derive this result does not provide an 
explicit decomposition of such states into product states, 
but merely proves that such a decomposition exists. It is. 



however, fairly simple to find such a decomposition |112] 
for mixed states with low values of e. 

The process begins by constructing an overcomplete 
basis for a single qubit; a basis set of this kind is sufficient 
to describe any state of a single qubit (a single spin) , but 
contains more basic elements than is strictly necessary. 
One suitable basis is the set of six states 
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which may be labeled as /3j, j = 1, 2, . . . 6. This can 
then be used to construct a two qubit basis set, (3jk, by 
taking direct products. This basis set, comprising 36 ele- 
ments, was constructed by taking products of single qubit 
states, and so is explicitly composed of product states 
only. Thus any density matrix which can be decomposed 
in this basis must be separable. 

As an example consider a pseudo-entangled state of 
the form given by equation p3 with e = 1/9: 
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This can be written in the form 






(58) 



and the proportions Pjk of each basis state Pjk can be 
obtained from the trace of the product of p and f3jk giving 
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This calculation also makes clear the danger of an un- 
questioning equation of entanglement and multiple quan- 
tum coherence. The state shown in equation p^ clearlyl 
contains double quantum coherence (it can be decom-' 
posed in product operator notation as \E/2-\-2IzSz/l-'&+\ 
DQx/^) and yet is not provably entangled. 
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